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ABSTRACT 

We report on Rossi X-ray Timing iSxpiorer observations of four type I X-ray bursters; namely 1E1724- 
3045, GS1826-238, SLX1735-269 and KS1731-260. The first three were in a low state, with 1-200 keV X- 
ray luminosities in the range ~ 0.05—0.1 LBdd (Lfidd: Eddington luminosity for a neutron star = 2.5x 10^* 
ergs s~^), whereas KS1731-260 was in a high state with a luminosity ~ 0.35 Lsdd- The low state sources 
have very similar power spectra, displaying high frequency noise up to ~ 200 Ifz. For KS1731-260, its 
power spectrum is dominated by noise at frequencies < 20 Hz; in addition a quasi-periodic oscillation at 
1200 Hz is detected in a segment of the observation. The 1-200 keV spectra of the low state sources are 
all consistent with resulting from thermal Comptonization with an electron temperature (kTs) around 
25-30 keV. For KS1731-260, the spectrum is also dominated by thermal Comptonization, but with a 
much lower kTs ~ 3keV and no significant hard X-ray emission. With the exception of GS1826-238, 
there is an underlying soft component, carrying at most ~ 25% of the total 1-200 keV luminosity. For all 
sources, we have detected an iron Ka line at 6.4 keV (although it is weak and marginal in 1E1724-3045). 
A reflection component is present in the spectra of GS1826-238 and SLX1735-269, and for both we find 
that the reflecting medium subtends only a small solid angle {Q./2-n ~ 0.15, 0.28). The origin of the line 
and the reflection component is most likely to be irradiation of the accretion disk by the X-ray source. 

We suggest a model in which the region of main energy release, where hard X-rays are produced would 
be an optically thin boundary layer merged with an Advection Dominated Accretion Flow (ADAF), 
and would be responsible for the rapid variability observed. The soft component observed probably 
represents the unscattered emission from an optically thick accretion disk of variable inner radius. When 
the accretion rate increases, the inner disk radius shrinks and the strength of the reflected component 
and associated iron line increase. At the same time, the Comptonization region cools off in response to 
an increased cooling flux from the accretion disk and from the reprocessed/reflected component, thus 
leading progressively to a quenching of the hard X-ray emission. If low state NSs accrete via ADAFs, 
the observation of X-ray bursts, indicating that all the accreting matter actually accumulates onto the 
NS surface, argues against the existence of strong winds from such accretion flows. 

Finally, we discuss two criteria recently proposed to distinguish between non-quiescent Black Holes 
(BHs) and Neutron Stars (NSs), and that are not contradicted by existing observations. The flrst 
one states that when thermal Comptonization is responsible for the hard X-ray emission, only BHs 
have kTE larger than ^ 50 keV. However, this criterion is weakened by the fact that there are NSs 
displaying non-attenuated power laws extending up to at least 200 keV, possibly implying non-thermal 
Comptonization or thermal Comptonization with kTE larger than 50 keV. The second criterion stipulates 
that only BHs are capable of emitting hard X-ray tails with 20-200 keV luminosities > 1.5 x lO^'^ergs 



1. Introduction 

It is now well established that hard X-ray emis- 
sion (E > 30 keV) from X-ray binaries is not exclu- 
sively associated with Black Hole systems (BHs). 
The major breakthrough came with the SIGMA 
and BATSE observations, which provided the first 
unambiguous detections of type I X-ray bursters 
(hence Neutron Star systems, NSs) at ~ 100 kcV 
(Barret & Vedrenne 1994, Barret & Tavani 1997 
and references therein). However, due to the mod- 
erate; sensitivity and low spectral resolution of 
these instruments and the lack of simultaneous X- 
ray observations, it was impossible to investigate 
the conditions under which NSs emit hard X-rays. 
Similarly, very little could be inferred about the 
accretion geometry and the relative contribution 
of the potential emitting regions to the total emis- 
sion (NS surface, boundary layer, accretion disk, 
corona). In addition, the hard X-ray data alone 
were not good enough to discriminate between the 
thermal and non-thermal models that have been 
put forward to account for this emission. Finally, 
if the emission of hard X-rays is indeed common 
to BHs and NSs, there may remain some differ- 
ences either in the shape of the hard tails (posi- 
tion of the energy cutoff; e.g. Tavani & Barret 
1997, Churazov et al. 1995), or in the luminosities 
that can be radiated simultaneously in soft and 
hard X-rays^ (Barret ot al. 1996, Van Paradijs 
& Van dcr Khs 1994, Barret & Vedrenne 1994, 
Zhang et al., 1996). All these potential differences 
have yet to be quantified precisely. They can now 
be addressed through observations performed by 
the Beppo-SAX and the Rossi X-ray Timing Ex- 
plorer (RXTE) satellites; the latter combining 
the Proportional Counter Array (PC A) and the 
High Energy X-ray Timing Experiment (HEXTE) 
(Bradt et al. 1993). These two instruments offer 
for the first time the possibility of observing NSs 
with good sensitivity and excellent timing capa- 
bilities simultaneously from 2 to ~ 150 keV. 

In this paper, we report on the RXTE ob- 
servations of four type I X-ray bursters; namely 
1E1724-3045, GS1826-238, SLX1735-269 and 
KS1731-260. All these sources have in common 
that they have already been detected up to ~ 100 

^In this paper, we define soft X-rays, as photons with energy 
between 1 and 20 keV (also X-rays), and hard X-rays, as 
photons with energy between 20 and 200 keV. 



keV at least once (Barret & Tavani 1997). The 
analysis has been made in a coherent way for all 
sources to allow a reliable and consistent compar- 
ison of their respective properties. The paper is 
organized as follows. First we review the rele- 
vant spectral observations of the four sources, and 
present their long term RXTE All Sky Monitor 
(ASM) light curves (Sect 2). We then describe 
our data reduction scheme, before presenting the 
results of our RXTE pointed observations (Sect 3). 
In section 4, we discuss the main features of our 
timing and spectral results and their implications 
for our understanding of NS accretion. Finally, we 
discuss the most recent observational criteria that 
have been proposed to distinguish non-quiescent 
BHs and NSs on the basis on their broad band 
spectral properties. 

2. The sources 

First, let us review briefly the previous obser- 
vations of these four sources. We will emphasize 
X-ray spectral observations, in particular previ- 
ous Nh measurements relevant to the PCA spec- 
tral analysis, existing hard X-ray observations, 
and finally more general information about these 
sources. 

2.1. 1E1724-3045 

1E1724-3045 is the most extensively studied of 
the four soiu'ccs discussed here. It is located in the 
globular cluster Terzan 2 in the general direction 
of the Galactic center. Its distance was first esti- 
mated to be 7.7 kpc (Ortolani et al. 1997), and 
recently revised to 6.6 kpc (Barbuy et al. 1998). 
It is a weakly variable X-ray source, and a per- 
sistent, though variable, hard X-ray source (Gold- 
wurm et al. 1993, 1994). Details about previous 
X-ray observations can be found in Barret et al. 
(1999a). Broad band spectral observations have 
already been performed by Beppo-SAX (Guainazzi 
et al. 1998). They showed that the 1 to 200 keV 
spectrum could be well fitted by the sum of a soft 
and a hard Comptonized component. The former 
could be equally fitted by a single temperature 
blackbody (BB) of fcrBB=0.6 keV, i?BB = 12 km, 
or a Multi-Color Disk blackbody (MCD Mit suda 
et al. 1984) with kTin=1.4 keV, Rin Vcos 6 ^=2.7 
km (at 10 kpc corresponding to Rin\/cos6'= 2.0 
km at 6.6 kpc). The hard Comptonized compo- 
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nent was fitted with a Com,ptt model (Titarchuk 
1994) with an electron temperature of ~ 27 keV, 
an optical depth ~ 3.3 for a spherical scattering 
cloud, and a temperature of 0.6 keV (kTvi/) for the 
seed photons. The ratio of the bolometric (0.1-100 
keV) luminosity of the soft component to the total 
luminosity was 11% and 36% for the BB and MCD 
models respectively. However, as pointed out by 
Guainazzi et al. (1998), whereas the radius de- 
rived from the BB fit is close to the NS radius, the 
very small value inferred for i?in^ appears unphys- 
ical (for any plausible values of the source inclina- 
tion); hence the hypothesis that the soft compo- 
nent comes from the boundary layer around the 
NS might be preferred. 

An ASCA observation confirmed the hardness 
of the source in X-rays, and allowed an accurate 
determination of Nh towards the source (~ 1.2 x 
10^^ H atoms cm^^). This value is consistent with 
that expected from the optical reddening of the 
cluster (Barret et al., 1999a). 

^Thc MCD model returns kTin and RinVcosS (for an as- 
sumed distance), which arc the color temperature of the in- 
ner accretion disk and the projected inner disk radius (note 
that in the case of a pseudo-newtonian disk with zero-stress 
inner boundary conditions around a Scliwarzschild BH, the 
actual inner disk radius is 2.73^^ times less than the mea- 
sured /tin, see Gierlinski et al. 1999). Correction for spec- 
tral hardening must be made to kTjn and Ri^ to account 
for the fact that the inner disk opacity is dominated by elec- 
tron scattering. We use a spectral hardening factor f=1.7, 
a value consistent with the source luminosity (Shimura & 
Takahara 1995), and determine the effective temperature 
as kTe£F=kTin/f, and the effective inner disk radius as 
Rgff= /2 _Rin= 2.9 Rin. For 1E1724-3045, this means that 
kTeff=0.82 keV, i?cff=11.4 km for an assumed inclination 
angle of 60° (the source is not a dipper, therefore its inclina- 
tion must be less than 75°). If we further assume that the 
disk terminates at the last stable orbit, then as described 
in Ebisawa et al. (1994), Rg (the Schwarzschild radius) is 
related to iJcff as 3R,s ~ 9 M>fs km oc ^ijR^fi (M^g mass 
of the neutron star in units of Mq , rj < 1 accounts for the 
decrease of RinVcos 9 by rclativistic effects, sec Ebisawa 
et al. 1994, here we take r] = 0.6 as in Shiimu a & Takahara 
(1995)). From this, the low inferred value of _R(,ff implies 
a very low mass of 0.4 Mq for the NS (for a 1.4 Mq NS, 
Rgff should be 37 km, equivalent to Rin ~ 13 km). This 
does not seem to favor the MCD model used to fit the soft 
component in the Beppo-SAX data of 1E1724-3045. How- 
ever, one has to keep in mind that ilin is derived from the 
normalization of the MCD model, and therefore would be 
underestimated if some fraction of the disk flux is scattered 
in a corona. The measured value of Hjn derived from the 
unscattered fraction would then depend on the geometry, 
and the optical depth of the scattering corona. 



2.2. GS1826-238 

When GS 1826-238 was discovered serendipi- 
tously by GINGA in 1988, it was first classified as 
a transient source, and further classified as a Black 
Hole Candidate (BHC) based on its rapid and in- 
tense flickering and hard Power Law (PL) spec- 
trum in X-rays (Tanaka 1989). The source was 
later detected by TTM (In't Zand 1992), ROSAT- 
PSPC (Barret et al. 1995a), and at hard X-ray 
energies by OSSE (Strickman et al. 1996). The 
BH nature of GS1826-238 was questioned (e.g. 
Strickman et al. 1996, Barret et al. 1996), but 
the firm evidence that GS1826-238 did not con- 
tain a BH came with the discovery of type I X-ray 
bursts (the unmistakable signature of a NS) with 
the Beppo-SAX Wild Field Gameras (WFC, Uber- 
tini et al. 1997). In the WFC data, the 70 bursts 
observed over more than 2.5 years of source mon- 
itoring exhibit a quasi-periodicity of 5.76 hours in 
their occurrence times (Ubertini et al. 1999). In 
addition, the so-called a parameter which is the 
ratio between the average persistent X-ray flux to 
the average flux emitted during X-ray bursts (as- 
suming that both are isotropic) was estimated to 
be 60 ± 7 for two of the 70 bursts analyzed by 
Ubertini et al. (1999). Similarly In't Zand et al. 
(1999), using also a sample of two bursts, observed 
a = 54 ± 5 consistent with the previous estimate. 
This value is consistent with a picture in which all 
the accreting material is accumulated onto the NS 
and helium is burned during the bursts (a should 
be ^ 20 and ^ 80 for pure Hydrogen and Helium 
burning respectively). 

As is the case for 1E1724-3045, broad band 
observations of GS1826 238 already exist. First, 
Strickman et al. (1996), combining non simul- 
taneous GINGA (X-ray) and OSSE (hard X-ray) 
data showed that the 2-200 keV spectrum of the 
source could be fitted by an exponentially cutoff 
power law (GPL: r= 1.76 ± 0.02, Ecutoff= 58 ± 5 
keV, r is the photon index of the power law, 
i.e. the energy index would be r-2) with a weak 
(and marginally significant) reflection component. 
More recently, Beppo-SAX observed the source, 
following a BATSE trigger which indicated that 
GS1826-238 was emitting hard X-rays (Frontera 
et al. 1998, Del Sordo et al. 1999). The 0.5- 
200 keV spectrum was well fitted by a compos- 
ite model consisting of a blackbody of 0.94 ± 0.05 
keV, and a PL of photon index T= 1.34 ± 0.04, 
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exponentially cutoff at 49 ± 3 keV. The column 
density measured by Beppo-SAX was 0.47 x 10^^ 
H atoms cm-^, consistent with the ROSAT PSPC 
value (0.5 ± 0.04 x 10^^ H atoms cm~^, Barret et 
al. 1995); both being slightly larger than the value 
expected from the optical reddening towards the 
source (E(B-V)^ 0.4 corresponding to 0.22 x 10^^ 
H atoms cm~^, Predehl & Smith 1995, Barret 
et al. 1995). Another Beppo-SAX observation 
performed half a year before by In't Zand et al. 
(1999) yielded best fit spectral parameters consis- 
tent with those reported in Del Sordo et al. (1999). 

GS1826 238 has a V= 19.1 ± 0.1 optical coim- 
terpart (Barret et al. 1995, Homer et al. 1998) 
in a possible 2.1 hour orbital period system, and 
is therefore a Low Mass X-ray Binary (LMXB). If 
the 2.1 hour modulation is the true orbital period 
of the system, one can try to estimate the source 
distance following the approach of Van Paradijs & 
McClintock (1994). In LMXBs, the bulk of the op- 
tical light originates from the reprocessing of the 
X-rays in the accretion disk. For systems for which 
reliable distance estimate exists, there is a good 
correlation between the absolute V magnitude and 
the X-ray luminosity {Lx) and the size of the ac- 
cretion disk (and hence the orbital period Porb, 
Lv oc L]['^ Porb^/^)- Most LMXBs have My in 
the range 0-2, whereas some short period systems 
have My in the range 3-5. With a 2.1 hr orbital 
period, GS1826-238 is a short period system, and 
if one assiimcs My in the above range (and an ab- 
solute dcrcddoncd V magnitude of 17.9), one gets 
a formal distance range between ^ 4 and ^ 10 kpc. 
Constraints on the source distance can also be in- 
ferred from the observation of X-ray bursts. In't 
Zand et al. (1999) derived an unabsorbed bolo- 
nictric peak flux of a burst of 2.7 ± 0.5 x 10~* ergs 
s^^cm^^. The fact that this burst did not show 
evidence for photospheric expansion implies that 
its luminosity is below the Eddington limit (LEdd) 
which we assume to be 2.5 x 10^* ergs s~^ (this 
value is appropriate for helium-rich material, a 1.4 
Mq NS, and a moderate gravitational redshift cor- 
rection. Van Paradijs & McClintock 1994). This 
in turn implies an upper limit on the source dis- 
tance of 9.6 kpc. In the remainder of this paper, 
we assume a distance of 7 kpc for GS1826-238. 



2.3. SLX1735-269 

SLX1735 269 was reported for the first time in 
1985, when it was detected by the Spacelab 2 X- 
ray Telescope (Skinner et al. 1987). However, it 
was present in the Einstein Slew Survey in ob- 
servations performed ~ 5 years earlier (Elvis et 
al. 1992). SLXl 735-269 was later detected by 
TTM (In't Zand 1992), by ART-P, ROSAT-PSPC 
(Grebenev et al. 1996) and ASCA (David et al. 
1997). SIGMA observations have shown that it 
is a persistent hard X-ray source of the Galactic 
center region (Goldwurm et al. 1996). The weak- 
ness of the source (~ 15.4 mCrab in the 35-75 
keV range) did not allow tight constraints to be 
put on the time averaged (1990-1994) hard X-ray 
spectrum, which could be fitted either by a simple 
PL of photon index 2.9 ± 0.3, or alternatively 
by a Comptonization model {Compst in XSPEC, 
Sunyaev and Titarchuk, 1980) with kTE of 26^^ 
keV (Goldwurm et al. 1996). Although early sus- 
pected to contain a NS based on the softness of 
its hard X-ray spectrum (Goldwurm et al. 1996), 
its nature remained uncertain until type I X-ray 
bursts were discovered with the Beppo-SAX WFC 
(Bazzano et al. 1997). The ASCA observations 
of SLX1735 269 revealed that its 0.6-10 kcV spec- 
trum could be well fitted with a PL of index 2.15 
(David et al. 1997), absorbed through an Nh of 
^ 1.4 — 1.5 X 10^^ H atoms cm~^, a value which is 
consistent with a source location near the Galactic 
center (i.e. at a distance of 8.5 kpc, David et al. 
1997). The Nh value derived with ASCA is also 
consistent with the one derived from the ROSAT- 
PSPC and ART-P observations (1.2 - 1.4 x 10^^ 
H atoms cm~^, Grebenev et al. 1996). The rapid 
variability of the source was recently investigated 
by Wijnands & Van der Klis (1999), using a short 
10 kilosecond Rossi X-ray Timing Explorer ob- 
servation performed between February and May, 
1997 (see below). 

No accurate distance estimate exists so far. 
From the burst reported in Bazzano et al. (1997a) 
(1.5 X 10~^ ergs s~^cm~^, 2-10 keV, Bazzano et 
al. 1997b), an upper limit of ~ 10 kpc can be de- 
rived (after bolometric corrections and assuming a 
blackbody of 2 keV for the burst). In the remain- 
der of the paper, we will assume a distance of 8.5 
kpc. 
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2.4. KS1731-260 

KS1731-260 was discovered in October 1988 by 
TTM (Sunyaev et al. 1990), and then classi- 
fied as a transient. KS1731-260 is located in the 
Galactic center region, about 5° away from the 
Galactic nucleus and only 1° from the X-ray pul- 
sar GXl-l-4. Above the mean persistent level of 
80 mCrab (2.7 x lO^i^ergs s-^cm^^ 1-20 kcV), 
TTM observed several X-ray bursts from KS1731- 
260, thus indicating that it contains a weakly mag- 
netized NS. The TTM spectrum was fitted by a 
Thermal Bremsstrahlung (hereafter TB'^) of 5.7 
keV, absorbed through an Nh of 2.2 x 10^2 H 
atoms cm^2. Follow-up ROSAT and optical ob- 
servations indicated that the source is a likely 
LMXB (Barret et al. 1998). The Nh derived from 
the ROSAT-PSPC all sky survey observations of 
1.3 X 10^2 H atoms cm"^ for a PL fit; is about a 
factor of 2 less that the value derived with TTM, 
thus possibly indicating intrinsic Nh variations 
within the source. More recently, we observed 
KS1731-260 with ASCA (September 27th, 1997) 
and confirmed the Nh found by ROSAT (Narita 
et al. 1999). Additional X-ray observations can 
be found in Yamauchi & Koyama (1990) and Alek- 
sandovich ct al. (1995). KS1731-260 was detected 
only once in hard X-rays between 35 and 150 keV 
by SIGMA, and has remained undetectable since 
then at these energies (Barret et al. 1992). Unfor- 
tunately, no simultaneous X-ray observations exist 
for the SIGMA detection. 

RXTE first observed KS1731-260 in July- 
August 1996. It was then in its standard high 
state. These observations led to the discovery of 
a highly coherent 524 Hz periodic X-ray signal at 
the end of the contraction phase of an X-ray burst 
that showed photospheric expansion (Smith et al. 
1997). This signal was tentatively interpreted as 
an indication of the NS spin (Smith et al. 1997). 
However, in the same data set, two simultane- 
ous High Frequency Quasi-Periodic Oscillations 
(HFQPOs) at 898 ± 3.3 Hz and 1158 ± 9 Hz were 

*Note that at such luminosities, TB fits are unphysical be- 
cause they imply emission measures 10®^ — 10^" cm^. The 
optically thin requirement leads to a size of 10-"' — lO-'-' 
cm for the cloud. These sizes are at least ~ 100 — 1000 
times larger than the region of main energy release around 
a NS. This means that one should be cautious about the 
Nj] fitted, because its value depends critically on the shape 
assumed to fit the continuum. 



also found (Wijnands & Van der Klis 1997), see 
below). The frequency separation (260.3 ± 9.6 Hz, 
corresponding to 3.8 msec) which is equal to half 
the frequency of the burst oscillations was then in- 
terpreted as the true NS spin frequency (Wijnands 
& Van der Klis 1997). 

In addition, assuming that the burst observed 
by RXTE reached the Eddington luminosity, a dis- 
tance of 8.8 ± 0.3 kpc was derived for the source 
(Smith et al. 1997). At this distance, the im- 
plied source X-ray luminosity at the time of the 
RXTE observations was 5.7 x 10^'^ergs s~^. Fi- 
nally a Nh of 6 X 10^^ H atoms cm"^ was derived 
for a TB fit of the RXTE PGA spectrum; a high 
value, significantly larger than the ASGA, ROSAT 
and TTM ones. 

2.5. ASM light curves 

The RXTE ASM long term light curves of 
all four sources are shown in Fig. 1 (for infor- 
mation about the RXTE/ ASM see Levine et al. 
1996). The mean ASM count rates are 2.9 cts 
s-i (- 38 mCrab), 2.9 cts s^^ (- 38 mCrab), 2.6 
cts s-i (~ 35 mCrab), 11.3 cts s'^ (~ 150 mCrab) 
in the 2-12 keV range for 1E1724-3045, GS1826- 
238, SLX1735-269 and KS1731-260 respectively. 
Assuming a PL spectrum of photon index 2, and 
the source Nh, at the source distance, these val- 
ues correspond to average 1-20 keV luminosities of 
~ 1.0 X 10^7 ergs 8'^= 0.04 L^dd, 8.8 x lO^^ ergs 
s-i= 0.03 LEdd, 1-3 X 10^'^ ergs 5-^= 0.05 Lsdd for 
1E1724 3045, GS1826 238 and SLX1735 269 re- 
spectively. For KS1731-260, assuming a thermal- 
like X-ray spectrum, the mean flux is about a fac- 
tor of 5 larger: ^ 4.9 x 10^"^ ergs s^^ = 0.2 LEdd- 
These light curves provide the first confirmation 
of the low variability of SLX1735-269, GS1826- 
238 and 1E1724 3045. In particular, for GS1826- 
238, this together with independent measurements 
(e.g. Barret et al. 1995; In't Zand et al. 1999) in- 
dicate that since its discovery in 1988, the source 
is a very stable accretor, and certainly not a clas- 
sical transient as previously thought. The light 
curve of KS 173 1-260, shows also that the source is 
persistent but occasionally undergoes low intensity 
states, which might be associated with episodes of 
hard X-ray emission. 
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Fig. 1.— The RXTE/ASM long term light curves 
of 1E1724-3045, GS1826-238, SLX1735-269 and 
KS1731-260. Time is given as days after January 
1st, 1994. Data have been selected using /select 
to satisfy the criterion that the measured rate for 
each individual "dwells" is positive and that the 
reduced of the fit to recover the source inten- 
sity is less than 1.5. Data from all three Scanning 
Shadow Cameras are combined within Icurve. The 
binning time of the light curves is 2 days. The gaps 
in the light curves are due to the sun getting too 
close to the source region (see Levine et al. 1996 
for information about the ASM) . The dates of our 
pointed observations are indicated with arrows. 

3. Observations and Results 

The PCA instrument consists of a set of 5 iden- 
tical Xenon proportional counters covering the 2- 
60 keV energy range with a total area of about 
6500 cm^ (Bradt et al. 1993). The HEXTE instru- 
ment is made of two clusters of 4 Nal(Tl) / CsI(Na) 
phoswich scintillation detectors providing a total 
effective area of 1600 cm^ in the 15-200 keV range''' 
(Rothschild et al. 1998). The two clusters rock 

*Note that this is about a factor of 4 larger that the effective 
area of the SIGMA imaging telescope for a source located 
in its fully coded field of view (Paul et al. 1991). 



alternately between the source and background 
fields to measure the background in real time. Spe- 
cial care has been taken to avoid the presence of 
any known X-ray sources in the background fields. 

The rapid variability of the 4 sources has been 
investigated with the PCA science event data. 
Our spectral analysis is based on the standard 2 
data for the PCA. They provide count spectra 
each 16 seconds in 128 energy channels covering 
the 2 to 100 keV range. The FTOOLS used for 
PCA background estimation is pcabackest version 
2.1b (released in October 1998). For HEXTE, we 
also use the standard mode data. They provide 
64 channel count spectra each 16 seconds. PCA 
spectra have been first accumulated for each PCU 
units for integration times varying from ~ 1000 
to ~ 2000 seconds, roughly equivalent to an or- 
bit after filtering for good time intervals. In par- 
ticular, PCA data recorded during and up to 30 
minutes after the SAA passage have been removed 
using the FTOOLS 4.2 version of xtefilt. PCA re- 
sponse matrices were then computed for each PCU 
using the FTOOLS pcarsp 2.36. We have then 
summed individual PCU spectra using addspec 
version 1.0.1, which combines not only the PHA 
files, but also the background PHA files, and com- 
putes the associated response matrices. 

For HEXTE, spectra were accumulated for each 
cluster with similar integration times to those used 
for the PCA. The latest response matrices have 
been used (hexte-97mar20c-pwa.rmf). Cluster A 
and B spectra were then combined together us- 
ing addspec, to get a single HEXTE spectrum av- 
eraged over the whole observation. The source 
and background light curves and spectra were cor- 
rected for dcadtime effects using version 0.0.1 of 
hxtdead. 

Version 10.00 of XSPEC (Arnaud 1996) has 
been used for the spectral fitting. When combin- 
ing PCA and HEXTE spectra, the energy ranges 
of the fit were 2.5 to 25 kcV for the PCA, and 
above 25 keV and up to at most 150 keV for 
HEXTE. These energy ranges were selected to ex- 
clude the steepening of the PCA spectrum above 
25 keV (in most spectra, it starts around 20 keV), 
as well as the flatening of the HEXTE spectrum 
towards low energies. Starting at 2.5 keV with the 
PCA implies that in most cases, the Nh cannot be 
very well constrained. In each case, we have there- 
fore set Nh to the most accurate value measured 
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so far (i.e by ASCA/ROSAT / Beppo-SAX). In or- 
der to account for the uncertainties in the relative 
calibration of the PCA and HEXTE experiments, 
we have also left the relative normalization of the 
two spectra as a free parameter of the model used. 

To investigate the systematics in the PCA data 
and to validate our analysis scheme, we have 
first analyzed a 15 kilo-second Crab observation 
(March 22, 1997). We have found that to get ac- 
ceptable fits, it was necessary to add a system- 
atic error to the data to reduce the effects of the 
imperfect knowledge of the instrument response 
near the K and L edges of xenon. These errors 
derived from looking at the residuals of the power 
law fit are 0.5% between 2.5 and 15 keV, 1% be- 
tween 15 and 25 keV and 2.5% above. Fitting 
the Crab spectrum between 2.5 and 25 keV thus 
yields a x^ ^f of 1.0 (50 d.o.f.), a power law in- 
dex of 2.18, and a normalization at 1 kcV of 13.5 
Photons cm-2 s'^ keV'^ (Nh= 3.3x lO^^H atoms 
cm~^ being frozen during the fit). The ratio be- 
tween the data and the power law folded through 
the response matrice is shown in Fig. 8. A sim- 
ilar level of systematics has been used by several 
groups (e.g. Rothschild ct al. 1999, Wilms et al. 
1999, Bloser et al. 1999). Based on this, before 
the fitting, we have combined quadratically to the 
Poisson errors of our data, this level of system- 
atics using the FTOOLS grrpha version 5.7.0. No 
systematic errors were added to the HEXTE data. 

3.1. The RXTE observations 

The RXTE observations are summarized in 
Fig. 2, 3, 4 and 5, for 1E1724-3045, GS1826- 
238, SLX1735 269, and KS1731-260. These fig- 
ures represent the HEXTE background-subtracted 
light curve (top panel), the PCA 2-40 background- 
subtracted light curve (underneath), the PCA 
color-color diagram (soft color: 7-15 keV/2-7 
keV, hard color 15-40 keV/7-15 keV (bottom left 
panel), and a PCA hardness intensity diagram 
(bottom right panel). For clarity and homogene- 
ity, the data from orbits contaminated by bursts 
have been filtered out (only data recorded with 
the 5 PCA units ON are shown). 

The observation of 1E1724-3045 started on 
November 4th, 1996 (TO^llh 54mn 39s), and was 
spread over more than 4 days. During the observa- 
tion, an X-ray burst occurred on Nov. 8th at 05h 
30mn 07s (UT). As said above for KS1731-260, a 



few NSs exhibit nearly coherent pulsations during 
X-ray bursts. Sampling the burst profile each sec- 
ond, we have searched for such a coherent signal 
in both the science event and burst catcher mode 
data. No significant signals were detected between 
200 and 1500 Hz. Data from the orbit when the 
burst occurred and the next one have been re- 
moved from the present analysis. The source did 
not display any significant variability over the ob- 
servation. In addition, looking at the color-color 
diagram, one can see that it was observed in a sin- 
gle spectral state (see also, Olive et al. 1998). A 
single PCA spectrum averaged over the whole ob- 
servation has thus been considered for the spectral 
analysis. 



lE17a4-3045 




0.530 0.540 0.550 0.560 0.570 
HRl (7-15 keV/a-7 keV) 



430 440 450 460 
COUNT RATE (2-15 keV) 



Fig. 2. - Top HEXTE 25-100 keV hard X-ray light 
curve with underneath the 2-40 keV PCA light 
curve of 1E1724 3045 over the 4.5 day span of our 
observation. Bottom left The color-color diagram 
(the soft color HRl is the ratio 7-15 keV/2-7 keV 
and the hard color HR2 is 15-40 kcV/7-15 keV) 
and Bottom right the hardness-intensity diagram 
for the soft color. Only PCA data recorded with 
the 5 PCA units ON are shown (this explains the 
first gap in the PCA light curve while HEXTE was 
recording data). 

GS 1826-238 was observed between November 
5th and 6th, 1997 for about 50 ksec also (T0=08h 
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22mn 39s, UT). There are two bursts in the middle 
of the observation. As for 1E1724-3045, no coher- 
ent pulsations were detected in those bursts. The 
time separation between the two bursts is ^ 5.6 
hours, consistent with the 5.76 hours (Ict spread 
of 0.26 h) periodicity reported by Ubertini et al. 
(1999). In addition, from this periodicity, one ex- 
pects a burst to have occurred ~ 1600 seconds be- 
fore the beginning of our observation. The tail of 
this burst is clearly seen in the first orbit of data. 
Data from this orbit, as well as from those contam- 
inated by the two bursts, have been removed from 
the present analysis. The spectral analysis of these 
bursts and their impact on the persistent emission 
will be reported elsewhere. As shown in Fig. 3 the 
source intensity varies by as much as ^ 10% in the 
2-40 keV range, but the source did not display sig- 
nificant spectral variability. As for 1E1724-3045, a 
single PCA spectrum was then used in the spectral 
analysis. In the hard X-ray band, it is the bright- 
est of the sources considered here (the count rate 
in HEXTE-Cluster A is 8.0 cts s'^ in the 25-100 
keV range). 

SLX1735-269 was observed in October 10th, 
1997 (T0=04h 37mn 19s); its observation which 
was scheduled for 50 ksec ended 2.5 days later. 
The mean PCA count rate is 165.5 cts s~^ (2-40 
keV), and this is the faintest of the four sources 
considered here. No X-ray bursts were observed. 
The source is clearly variable on minute time 
scales. However it is clear that these intensity 
variations are not accompanied by strong spectral 
changes, as the source remains in a limited region 
in the color-color diagram. 

Finally, the observation of KS1731-260 started 
on October 28th, 1997 (T0=22h 20mn 15s, UT), 
lasted for about 50 ksec, and ended on October 
30th. In the PCA, KS1731-260 has the largest 
count rate of the sources considered here. The so- 
called banana shape is clearly visible on the color- 
color and hardness intensity diagrams. It moves 
from the lower to the upper parts of the banana 
as the count rate increases (by up to ~ 25%, 2-40 
keV) . Given the source variability, we have consid- 
ered three PCA spectra in the spectral analysis; 
one for each day of the observation. It is detected 
by HEXTE only up to ~ 50 keV. No X-ray bursts 
were detected. 



GS1626-23B 
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Fig. 3.— Top HEXTE 25-100 keV hard X-ray light 
curve with underneath the 2-40 keV PCA light 
curve of GS1826-238 over the 1.2 day span of our 
observation. Bottom left The color-color diagram 
(the soft color HRl is the ratio 7-15 keV/2-7 keV 
and the hard color HR2 is 15-40 keV/7-15 keV) 
and Bottom right the hardness-intensity diagram 
for the soft color. The data contaminated by X- 
ray bursts have been removed. Only PCA data 
recorded with the 5 PCA units ON are shown. 

3.2. Timing properties 

The normalized Power Density Spectra (PDS) 
averaged over the whole observation of all four 
sources are shown in Fig. 6. The first three sources 
show noise up to ~ 200 Hz. The integrated power 
of these PDS is 29.1%, 26.1%, 27.6% in the 2- 
40 keV band (0.005-300 Hz). In addition, they 
also show a break in the PDS at low frequen- 
cies (typically around 0.1-0.2 Hz), and a QPO- 
like feature around 1 Hz. Although the complete 
analysis is still under way (for GS1826-238 and 
SLX1735-269), no strong HFQPOs were detected 
from any of these three systems. For 1E1724 3045, 
for which a sensitive search has already been con- 
ducted, an upper limit of 2.5% on the RMS of a 
1000 Hz QPO has been derived (5-30 keV) (Barret 
et al. 1999b). 
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Fig. 4.— Top HEXTE 25-100 keV hard X-ray light 
curve with underneath the 2-40 keV PCA light 
curve of SLX1735-269 over the 2.5 day span of our 
observation. Bottom left The color-color diagram 
(the soft color HRl is the ratio 7-15 keV/2-7 keV 
and the hard color HR2 is 15-40 keV/7-15 keV) 
and Bottom right the hardness-intensity diagram 
for the soft color. Only PCA data recorded with 
the 5 PCA units ON are shown. 



Fig. 5.— Top HEXTE 25-50 kcV hard X-ray light 
curve with underneath the 2-40 keV PCA light 
curve of KS1731-260 over the 1.5 day span of our 
observation. Bottom left The color-color diagram 
(the soft color HRl is the ratio 7-15 keV/2-7 keV 
and the hard color HR2 is 15-40 keV/7-15 keV) 
and Bottom right the hardness-intensity diagram 
for the soft color. Only PCA data recorded with 
the 5 PCA units ON are shown. 



The PDS shown in Fig. 6 for SLX1735 269 is 
very similar both in shape and normalization to 
the one reported by Wijnands & Van der Klis 
(1999a) from a short RXTE observation per- 
formed in February-May 1997. Although the 
source count rate did not change a lot during our 
observation (see Fig. 4), we have computed two 
PDS for segments of the observation associated 
with the highest and lowest count rates: 172 cts 
s~'^ (156 cts s"'^) versus 162 cts (150 cts s^-^), 
2-40 keV (2-16 keV), respectively. Fitting these 
two PDS with broken power laws yields a break 
frequency of 0.15 ± 0.1 Hz, and 0.08 ± 0.02 Hz 
respectively. Thus, within our limited range of in- 
tensity variations, the break frequency decreased 
with the count rate, a behaviour which is generally 
observed from atoll sources (Van der Klis, 1994). 
At lowest count rate, Wijnands & Van der Klis 
(1999a) reported the opposite behaviour between 



two segments of observations of count rates 141 
and 113 cts s~^ (2-16 keV); i^Break increased from 
0.11 to 2.3 Hz. Therefore, this means that like 
in the case of the millisecond pulsar SAXJ1808.4- 
3658 (Wijnands & Van der Klis 1998), it appears 
that ^'Break Correlates with the source count rate 
(the latter most likely tracking the mass accre- 
tion rate), down to a certain value, below which 
I'Break and the source count rate anti-correlates. 
For both SLX1735 269 and the millisecond pul- 
sar, this transition occurred at a luminosity of 
~ 2 - 4 X 103<5 ergs s'^ (3-25 keV). 

KS1731-260 is characterized by weaker variabil- 
ity with an integrated power of ~ 3.0% (2-40 keV). 
Its PDS can be approximated by a PL of index 
a = 1.3 on which a QPO centered at 7.6 ± 0.7 Hz 
is superposed (FWHM = 3.1 ± 2.0 Hz for a gaus- 
sian fit). This QPO is very similar to the so-called 
"Horizontal branch" QPO observed in Z sources. 
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Fig. 6. — The power density spectra of all four 
sources averaged over the whole observations. 
1E1724-3045, GS1826-238 and SLX1735-269 dis- 
play high frequency noise extending up to 200-300 
Hz with RMS amplitudes of 29.1%, 26.1%, 27.6% 
in the 2-40 kcV band (0.005-300 Hz). KS1731- 
260 displays Very Low Frequency Noise with a 
RMS amplitude of ~ 3.0% (2-40 keV). These PDS 
are consistent with the first three sources being in 
the so-called island state (low state), and KS1731- 
260 being in the banana state (high state). 

In addition, at higher frequencies, for the segment 
of the observation performed on October 28th (i.e. 
at the lowest count rate ^ 1050 cts s^^, 5-30 keV), 
we have detected at the 4.7cr level a HFQPO cen- 
tered at - 1200 ± 10 Hz (FWHM=45.0 ± 20 Hz, 
RMS=2.5 ± 0.4%)5. This is the highest HFQPO 
ever detected from KS1731-260. No HFQPOs were 
detected in the subsequent observations performed 
on October 29th (1285 cts s^^) and 30th (1355 cts 
s^^). We have derived an upper limit of ^ 2% on 
the RMS of any HFQPOs around 1000 Hz for the 

^The 1200 Hz HFQPO observed is a lower limit of the Keple- 
rian frequency at the innermost stable circular orbit (Kluz- 
niak 1998, Miller et al. 1998). This translates into an upper 
limit on the mass of the NS: Mns < 1-83 X {l + 0.75j) M©, 
where j = cJ/GM'^ is a dimensionless stellar angular mo- 
mentum, and 0.1 < j < 0.3 for a NS spin frequency of 
~ 260 Hz (KluzniaJj 1998). This yields an upper limit of 
2.25 Mq for the NS in KS1731-260. 



October 29th and 30th observations, respectively. 
Our October 28th detection follows the correlation 
between the frequency of the HFQPO (i'hfqpo) 
and the count rate, as well as the anticorrelation 
between fHFQPO and RMS presented in Wijnands 
& Van der Klis (1997). 

3.3. Spectral fitting 

Previous observations of similar NSs have 
shown that their broad band spectra can be gen- 
erally fitted by the sum of two components (Mit- 
suda et al. 1984, White et al. 1988, Mitsuda et 
al. 1989). The first one is soft, contributes to 
the spectrum mainly below 10 keV, is modeled 
by a multi-color disk blackbody or a single tem- 
perature blackbody, and is interpreted as coming 
either from the accretion disk, or the NS surface, 
or an optically thick boundary layer (e.g. White et 
al. 1988, Mitsuda et al. 1984). The second one is 
harder and dominates the spectrum above 10 keV; 
it is often modeled by simple PLs (below 20 keV). 
However, when high energy coverage exists above 
30 keV (e.g. with Beppo-SAX and RXTE), comp- 
tonization models provide more physical fits (e.g. 
Guainazzi et al. 1998, In't Zand et al. 1999). The 
comptonization process is speculated to take place 
in a scattering corona located somewhere in the 
system; around the NS (e.g. optically thin bound- 
ary layer, spherical corona) or above the disk. In 
some cases, a relatively strong iron Ka line (6.4 
keV) is also observed above the continuum (e.g. 
White et al. 1986). In a few cases (e.g. 4U1608- 
522, Yoshida et al. 1993, SAXJ1808.4-3658, Gil- 
fanov et al. 1998), such a line is accompanied by a 
broad absorption like feature, interpreted as par- 
tial absorption or refiection by a cold or weakly 
ionised medium of the intrinsic PL component. 
Hereafter we report on the results of the spectral 
analysis of our observations, placed in the frame- 
work of these previous results. 

3.3.1. 1E1724-3045 

Preliminary results of the analysis of the 
1E1724 3045 spectral data can be found in Olive 
et al. (1999). In the present analysis, Nh has been 
set to the ASC A/ Beppo-SAX value (i.e. 1.2 x 10^^ 
H atoms cm~^). Driven by the recent Beppo- 
SAX results, we have found that the broad band 
spectrum of 1E1724-3045 can be adequately de- 
scribed by the sum of two components; a soft 
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1E1724 3045 



Parameter 


C+BB+L 


C+MCD+L 


T 

fcTBB, kTin 
-RbB, RinVcOS^ 

EqW 


00 1+3.0 

28-113.0 

-"^•-■^-o.i 

lO.l+i'? 
2lt^° ■ 


25.6+^:i 

1 c+0.2 

J^-"-0.2 
q q+0.3 
"^•^-0.3 

'-■■^-0.2 
0+0.2 

27^1° 


(d.o.f) 


79.9 (80) 


85.7 (80) 


Fl_20 keV 


1.52 


1.58 


F2O-2OO keV 


0.90 


0.91 


/bb (%) 


14.4 


27.2 



Table 1: Best fit spectral results for 1E1724-3045. Nh has been set to 1.2 x 10^^ H atoms cm-^. The 
C+BB+L model is Com;)tt+Blackbody+a 6.4 keV narrow line (fTFe=0.1 keV), whereas C+MCD+L is 
Comptt+Multicolor Disk blackbody+a 6.4 kcV narrow line. kTE is the electron temperature in kcV. kTw is 
the temperature of the seed photons in the Comptt model (in keV). r is the optical depth of the spherical 
scattering cloud. ^Tbe is the temperature of the blackbody, kTjn is the inner disk color temperature derived 
from the MCD model. i?BB is the equivalent radius of the blackbody, whereas Rin\/cos^ is the projected 
inner disk radius (Mitsuda et al. 1984); they are both given in kilometers and scaled at the source distance 
(6.6 kpc). EqW is the equivalent width of the iron line in eV. The 1-20 keV flux (Fi_2o kev) and hard 
X-ray 20-200 keV flux (F20-200 kov) arc given in imits of xlO^^crgs s^-'^cm"^. /bb is the contribution of 
soft component (BB or MCD) to the total luminosity computed in the 1-200 keV range. Errors in all the 
tables are quoted at the 90% confidence level (x^=X^+2-7). 



and a hard Comptonized component well fitted 
by the Comptt model in XSPEC (Titarchuk 1994, 
Guainazzi et al. 1998, see also Barret ct al. 1999). 
The soft component is well described by a BB 
(X^=88.1, 81 d.o.f). A MCD could fit it as well, 
although the is slightly larger (x^=92.3. 81 
d.o.f). For the Comptonizing cloud, we have de- 
rived an electron temperature (kTs) of ~ 27 — 30 
keV and an optical depth of ^ 2.8 — 2.9 and ^1.1 
assuming that it has a spherical or disk-like ge- 
ometry respectively. In both cases, a temperature 
of ^ 1.0 kcV for the seed photons was derived. 
These parameters, which are listed in Table 1, are 
strikingly close to those found by Guainazzi et 
al. (1998), suggesting that in addition to its low 
variability in intensity (see Fig. 1), the source is 
also stable from the point of view of its energy 
spectrum. For the soft component, the best fit 
parameters are strikingly similar to those derived 
from the Beppo-SAX observations. 

Independently of the model fitting the contin- 
uum, an excess around 6 keV seems to remain 
in the residuals. Assuming that it results from 
iron K fluorescence (6.4 keV), and that it is nar- 



row (i7Fe=0.1 keV), the Comptt+BB+Line model 
yields a of 79.9 (80 d.o.f). The line equivalent 
width is 21 eV. In the remainder of this paper, 
we will assess the significance of additional com- 
ponents in the fitted model by means of the stan- 
dard F-test. What we call the F-test probability 
will always refer to the probability of rejecting the 
correct hypothesis, which states that, with the ad- 
ditional component the fit is better. We define 
Fm as {xl - xDli^i - i^2)/(xi/i^2) with xl > xl 
and v\ > V2 (where Xi 2 are the Chi-squared val- 
ues of the two fits, with t^i.2 number of degrees of 
freedom, d.o.f.); Fm follows a Snedecor's F distri- 
bution. In the present case, we have Fm=8.2, and 
P{F > _F'„i)=0.005, which means that the addition 
of the line to the Comptt+BB model is significant 
at the level of ~ 99.5%. It is worth pointing out 
that the significance of any features in the spec- 
trum depends strongly on the level of systematics 
added to the data, and that there are some uncer- 
tainties in the estimate of these systematics. To 
illustrate this, if instead of using 0.5% around 6 
keV, one uses 1%, as for example in Rothschild et 
al. (1999), the significance of the line would drop 
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Fig. 7. — The data and folded model (top panel), 
the residuals (center panel) and the unfolded spec- 
trum (bottom panel, Comptonized (Comptt) com- 
ponent is the dashed line, multi-color disk black- 
body model is the dot-dashed line, the Gaussian 
Une is the dotted Une) of 1E1724-3045. 

down to 85%, and hence the line detection would 
not be significant. Leaving the line energy and its 
width as free parameters does not improve the fit 
significantly. The presence of the line is also signif- 
icant at the level of 98.5% for the Comptt+MCB 
model. However, since no such a line was found in 
ASCA and Beppo-SAX observations, wo consider 
our detection very marginal. Fig. 7 shows the data 
and folded model {Comptt+BB+L) , the residuals 
and the unfolded spectrum of 1E1724-3045. 

For comparison with SIGMA data (35-200 
keV), we have fitted the HEXTE data alone in 
the 25-150 keV range. Although the fit is not sat- 
isfactory due to the presence of a clear high energy 
cutoff (~ 70 keV), a PL fit yields a photon index of 
2.7 ±0.1, to be compared with the time-averaged 
value observed by SIGMA (3.0±0.3, Goldwurm et 
al. 1993, 1994). Fitting the 2.5-25 keV PGA spec- 
trum with a simple PL (note an acceptable fit) 
yields a photon index of ^ 2.0. Wc therefore con- 
clude that the steepness of the hard tail observed 



by SIGMA was artificial, and due to the pres- 
ence of a high energy cutoff around 60-70 keV. We 
have also fitted the continuum with the relativis- 
tic Comptonization model developed by Poutanen 
& Svensson (1996) [CompPS in XSPEG). For the 
Comptonized tail, this }delds best fit parameters of 
35 keV for kTs and 2.1 for r (spherical geometry 
assumed). 

Prom this observation, we derive Li_2o kev= 
8.1 X 10^6 ergs s^^, and L20-200 keV= 4.8 x lO^^ 
ergs s~^ (d=6.6 kpc). Li_2o keV is consistent with 
the time-averaged value derived from the ASM 
light curve. 

3.3.2. GS 1826-23 8 

For GS1826-238, we have again set Nh to the 
value observed by ROSAT and Beppo-SAX (0.5 x 
10^^ H atoms cm~^). To illustrate the complexity 
of the spectral shape, in Fig. 8, we show the ra- 
tio between the PGA data and a PL model folded 
through the PGA response matrix. There is a clear 
excess around 6.0 keV. First, of the single compo- 
nent models, the GPL is the one that provides 
the best fit. This yields a cutoff energy at ~ 95 
keV and a photon index of 1.7 (x^=195.0 for 87 
d.o.f). A fit using a comptonization model (e.g. 
Compst in XSPEG) yields kTg of 20 keV and an 
optical depth of 4.6. Adding a narrow 6.4 keV line 
(o"Fc=0.1 keV) to account for the feature of Fig. 8 
improves the fit (x^=149.4 for 86 d.o.f) but still, 
the residuals show systematic deviations that in- 
dicate that this model is not the right description 
of our data. 

Since a reflected component had been found in a 
combined analysis of the OSSE and GINGA data 
(Strickman et al. 1996), we have tested for the 
presence of such a component by adding an ab- 
sorption edge at 7.1 keV. This leads to a significant 
decrease of (x^=100.5 for 85 d.o.f, Fm=Al.2>, 
P{F > Fra)> 7.9 X 10-9). Thus we have substi- 
tuted the GPL model by the XSPEG Pexrav model 
(Magdziarz and Zdziarski 1995), which is the sum 
of a GPL plus a Compton reflected component. 
The inclination angle to the source is unknown, 
but since GS1826-238 is not a dipper, we have as- 
sumed a standard value of ^ = 60°, and we leave 
the reflection scaling factor as a free parameter (we 
assumed solar abundance for the reflecting ma- 
terial). This yields a reflection scaling factor of 

0.20, and a narrow 6.4 keV line (crFe=0.1 keV) 
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GS1826-238 


Parameter 


CPL+L 


CPL+R+L 


r 

Ecutoff 
frefl 

Line Energy 

<7Fe 

EqW 
(d.o.f) 


1 7Q+^-^^ 
^■1^-0.02 

6.4 (fixed) 
0.1 (fixed) 

149.4 (86) 


1 79+0.01 

J--'^_0.01 

90.0+?[! 

"■^^-0.04 

6.11^2 

0.48l°:g 

85.2 (83) 


Fi_20 keV 


1.43 


1.44 


F2O-2OO keV 


1.17 


1.16 



Table 2: Best fit spectral results for GS1826 -238. Nr has been set to 0.5 x 10^^ H atoms cm'^. The CPL+L 
model is the sum of a Cutoff Power Law (P is the photon index, Ecutoff is the cutoff energy in keV), plus 
a narrow 6.4 keV line (crFe=0.1 keV). The CPL+L+R is the same but with the reflected component added 
(the CPL is substituted by the Pcxrav model in XSPEC, Magdziarz and Zdziarski 1995), and the line energy 
and width (crpe) left as free parameters of the fit. frefl is the reflection scaling factor normalized to the CPL. 
It should be equal to 1 for an isotropic source above an infinite flat disk. Same units for Fi_2o keV and 
r 20-200 keV aS in Table 1. 



of EqW of 37 eV (x^=97.1, 85 d.o.f.). The sig- 
nificance of the refiection component is very high 
(i^„=45.4, P{F > Fm)> 1.8 X 10-9). Leaving 
the line energy and its width as free parameters 
of the model, we get x2=85.2 (83 d.o.f.) with the 
line parameters: crFe=0.48 keV, centroid energy 
at &.ltl\ keV (90% confidence level, stiU con- 
sistent with a fluorescent iron Ka line), and an 
EqW of 50 eV. This decrease of is significant at 
the level of 99.7% (i^m=6.0, P{F > j;„)=0.003), 
we therefore conclude that this is the best fit to 
our data. As a further step, we have tried a re- 
flection model including ionisation: Pexriv model 
in XSPEC (Magdziarz and Zdziarski 1995). This 
model fits the data as well as the Pexrav model 
(x^=95.5, 84 d.o.f.), and yields a disk ionisation 
parameter consistent with 0. Therefore, with both 
models, our data are consistent with reflection 
from a cool neutral medium. 

We have also reanalyzed the GINGA spectrum 
of GS1826-238 used in Strickman et al. (1996) and 
Zdziarski et al. (1999). The reflection component 
is highly significant, and consistent with coming 
also from a neutral medium also. A fit with a sim- 
ple power law and a 6.4 keV narrow line yields 
X^=39.6 for 31 d.o.f., while adding the reflection 
component yields x^=16.1 for 29 d.o.f. {Fm=21.2, 
P{F > fJ)=2.1 X 10"^). The best flt parame- 
ters are r=1.88 ± 0.05, and frefl=0.8 ± 0.3, and an 



equivalent width for the iron line of 15 eV^. 

There is no soft component in the source spec- 
trum. With the Pexrav model, we have set a 90% 
confidence limit of about 1% on the fraction of 1- 
200 keV luminosity in the soft component modeled 
by a 1.5 keV BB component. The results of the 
best fit using the CPL-|-Line-|- Reflection model are 
given in Table 2. The data and folded model, the 
residuals and the unfolded spectrum of GS1826- 
238 are shown in Fig. 9. 

For comparison with hard X-ray observations of 
similar systems, we have fitted the HEXTE spec- 
trum with a PL. Although the fit is poor due to 
the presence of a clear cutoff in the spectrum, one 
gets a photon index of ^ 2.3 ± 0.2, similar to the 
value derived for 1E1724-3045. Fitting the hard 
tail with the CompPS model yields kTE=41 keV 
and r = 1.9. 

For this observation, one derives Li_2o keV= 



The evidence for a reflection component is very strong in 
both the RXTE and GINGA data. It is therefore quite 
puzzling that no such component has been reported so far 
from the two Beppo-SAX observations performed, while 
GS1826-238 was at a similar intensity level (7.7 x IQ-io 
ergs s" '^cm— ^ for RXTE, as opposed to 5.5 X 10" ^" ergs 
s-lcm-2 for Beppo-SAX (2-10 keV), Del Sordo ct al. 1999, 
In't Zand et al. 1999). Wc arc currently investigating this 
issue, searching for instrumental effects, sensitivity limits, 
and systematic errors in the fitting the continuum shape. 
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(d.o.f) 


45.2 (59) 


28.4 (57) 


Fl_20 keV 


0.72 


0.75 


F2O-2OO koV 


0.39 


0.30 


/bb(%) 


9.0 


12.6 



Table 3: Best fit spectral results for SLX1735-269. Nh has been set to 1.5 x 10^^ H atoms cm-^. The first 
model is the sum of CPL+MCD and a 6.4 keV narrow line. The best fit model is the sum of a CPL and 
its Compton reflected component (R) allowing for ionisation (Pcxriv in XSPEC), a MCD, and a narrow 6.4 
keV line (L) {ape=0.1 keV). The cutoff energy in Pexriv has been set to 200 keV and frozen during the 
fit. is the ionisation parameter in Pexriv. The contribution for the Galactic diffuse component has been 
modeled by a Raymond-Smith model of 3.4 kcV (normal abundance) and normalization factor of 1.5 x 10~^ 
in XSPEC units. See Table 1 for definition of /bb and units of Fi_2o keV and F20-200 keV. 



8.4 X 10^6 ergs s-\ and L20-200 kcV= 6.9 x IQ-''^ 
ergs s~^ (d=7 kpc). Li_2o keV is very close to the 
time-averaged value derived from the ASM light 
curve. 

3.3.3. SLX 173 5-269 

A preliminary report of the spectral analy- 
sis of SLX1735-269 can be foimd in Skinner et 
al. (1999). As in the cases of GS1826-238 and 
1E1724-3045, SLX1735-269 was observed in a 
hard state. However, it was a factor of two fainter 
than the other two sources, so that the statis- 
tical quality of the HEXTE spectrum is not as 
good. For that reason, we considered the HEXTE 
data only up to 50 keV. In our fits we fix Nr at 

1.5 X 10^^ H atoms cm~^, found in David et 
al. (1997), as the ASCA data are more sensi- 
tive to this parameter. We have first fitted the 
PCA-I-HEXTE spectrum of SLX1735-269 with a 
simple PL. The fit is not good (x^ ^ j > 6.5) but 
shows clear evidence for a broad feature around 6.4 
keV (see Fig. 8). As no iron line was seen in the 
ASCA data, and because of the relative faintncss 
of the source, we have first considered the possi- 
bility that the line might not be from the source 
but might represent the flux within the PC A 1° 
field of view from the galactic bulge diffuse emis- 



sion. We have therefore included in our model a 
contribution from the spectrum given by Valinia 
& Marshafl (1998) for their "Region 2" which in- 
cludes this part of the sky. Refltting the data with 
this component affects the shape of the residuals at 
low energies (< 5 keV); they now clearly show the 
presence of a soft component and reveal that the 
broad feature between 6-7 keV cannot all be ac- 
counted for by the diffuse emission, meaning that 
there remains a need for a signiflcant contribution 
from the source to the line. Such a line was also 
needed in the data analyzed by Wijnands & Van 
der Klis (1999). The soft component was fitted 
with a MCD (although the BB model fits the data 
equally well), whereas the line contribution from 
SLX1735 269 was modeled by a gaussian centered 
at 6.4 keV. This leads to a of 45.2^ (59 d.o.f), 
and thus the addition of the line and soft com- 



'''For SLX1735-269, the associated with the best fits are 
lower than 1 (sec Table 3). This makes naturally question- 
able the results of our F-test. Similarly, the errors com- 
puted on the best fit parameters should not be considered 
as true 90% uncertainties (these errors are computed un- 
der the assumption that the errors on the data are Poisso- 
nian). The low indicates certainly that for SLX1735— 
269, which is the faintest of the 4 sources, the systematics 
assumed are too large. However, for consistency in the 
analysis, we have chosen to set them to the values used for 
the other sources. 
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Fig. 8. — The ratio between the data and a folded 
power law model for the Crab, SLX1735-269 and 
GS1826-238. The response matrice used in the 
Crab analysis has been obtained using the same 
procedure as for the other sources. For GS1826- 
238 and SLX1735~269, this plot indicates the pres- 
ence of a strong feature around 6-7 keV. 



ponent is statistically significant based on our F- 
test. The EqW found (67 eV if the line is narrow, 
crFe=0.1 keV) is well below the upper limits of 150 
eV (c7Fe=0.1 keV) derived by David et al. (1997) 
using ASCA SIS data. 

However, by looking at the residuals above ^ 7 
keV the data show systematic deviations with an 
edge-like shape. An edge around these energies is 
suggestive of the presence of a reflection compo- 
nent. We have therefore first substituted the PL 
model with the Pexrav model in XSPEC (assum- 
ing no energy cutoff in the model, and a 6.4 kcV 
line with a crFe=0.1 keV, an inclination angle of 
60°). This leads to x^=31-4 (58 d.o.f.), indicating 
that the presence of a reflected component is sig- 
nificant at a level greater than 99.99% (^"^=23.1, 
P{F > Fm)= 1.1 10"^). The solid angle inferred 
is ^ 0.26 and the line EqW is 47 cV. Whether or 
not the reflector is ionised has been tested by us- 




10 100 
Energy (keV) 

Fig. 9. — The data and folded model (top panel), 
the residuals (center panel) and the unfolded spec- 
trum (bottom panel) of GS 1826-238. The iron line 
at 6.4 keV is represented by the dashed line, the 
reflection component by the long dashed line. 



ing the Pexriv model. The statistical quality of the 
data does not allow us to determine the line energy 
simultaneously with the Compton reflected com- 
ponent. We have assumed the line at 6.5 keV (ex- 
pected from a moderately ionised medium). This 
yields a x^=28.4 (57 d.o.f.), an ionisation param- 
eter, albeit poorly constrained, of ~ 7017q", a 
solid angle of ~ 0.28 and a line EqW of 39 eV 
(see Table 3). The improvement of from the 
Pexrav to the Pexriv model is significant at the 
level of 98.3% (F„=8.8, P{F > Fm)= 4 x lO'^). 
It is worth pointing out that leaving out the line 
at 6.4 keV increases the significance of the ionisa- 
tion of the reflector, but then the model becomes 
unrealistic, as the line is not expected to be at 6.4 
keV for the ionisation parameters derived from the 
flt (^, > 70). To conclude, one can say that our 
data for SLXl 735-269 are consistent with reflec- 
tion from a neutral or moderately ionised medium. 

In the absence of any cutoffs in the useful en- 
ergy range of the present measurements, comp- 
tonization models cannot really be used. We note 
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Fig. 10.— The data and folded model (top panel), 
the residuals (center panel) and the unfolded spec- 
trum (bottom panel) of SLX1735-269. The power 
law component is represented by the dashed line, 
the multi-color disk blackbody model by the dot- 
dashed line, the iron 6.4 keV line by the dotted line 
and the reflection component by the long-dashed 
line. The Raymond-Smith thermal spectrum cor- 
responding to the galactic diffuse emission around 
SLX1735-269 is also indicated as a dot-dot-dot 
dashed line. 

however that a Comptonization model for which 
the electron temperature would be ^ 30 keV (as 
seen in 1E1724-3045) could fit the continuum as 
well as the PL. This would yield an optical depth 
of ~ 3 for the scattering cloud. Similar results 
were obtained by combining simultaneous ASCA 
and SIGMA data (Goldwurm et al. 1995). The 
best fit parameters arc listed in Table 3 while the 
data and folded model, the residuals and the un- 
folded spectrum are shown in Fig. 10. 

For this observation, we derive Li_2o koV= 
6.0 X 10^6 ergs s-\ and L20-200 keV= 3.5 x lO^^ 
ergs s~^ (d=8.5 kpc). Li_2o keV is about a factor 
of two lower than the time-averaged value derived 
from the ASM light curve. 



3.34. KS 1731-260 

KS1731-260 was observed in a high state, and 
its spectrum is clearly much softer than the other 
three sources examined before. We have assumed 
the Nh measured by ROSAT and ASCA (i.e. 
1.2 X 10^^ H atoms cm~^). As the source inten- 
sity increased smoothly during our 3-day observa- 
tion, and its intensity variations are accompanied 
with spectral variations (see Fig. 5), as said above 
we have made a PCA spectrum for each of the 
three days. These PCA spectra were fitted simul- 
taneously with a single HEXTE spectrum aver- 
aged over the whole observation to increase the 
statistics at high energies. 

We first analyzed the October 28th spectrum. 
As there is a clear cutoff in the spectrum around 
10 keV, we have first fitted the continuum with 
a simple comptonization model (Compst). This 
model alone is clearly rejected {x^> 977 for 55 
d.o.f). This is mainly due to the presence of a 
broad feature again centered around 6.4 keV. We 
have therefore added an iron K fluorescence line 
(6.4 keV) and fitted its width and intensity. This 
improves the fit significantly (x^= 156.2 for 53 
d.o.f). However, looking at the residuals indicates 
that the low energy part of the spectrum is not well 
accounted for. As in similar systems, a soft com- 
ponent (either a blackbody or a disk blackbody) is 
usually observed (e.g. White et al. 1988), we have 
included such a component in the fit. This again 
improves the fit (x^=48.6 for 51 d.o.f., Fm— 56.5, 
P{F > Fm)=l.2 X 10-13) foj. addition of the BB 
model. We have tried to fit the "hard" compo- 
nent with a MOD model (Mitsuda et al. 1984), 
instead of using Compst, but the fit is rejected, as 
the reduced exceeds 4. On the other hand, the 
soft component could equally well be fitted with 
a MCD. For both models, the equivalent radius 
of the BB and the projected inner disk radius of 
the MCD are small; Rbb^ 5 km and RinVcos 9< 
4 km respectively, smaller than the NS radius, or 
the expected inner disk radius around a NS (espe- 
cially if the NS is as massive as 2Mq, see above). 
This does not mean however that the models can 
be ruled out, because for instance only a small 
fraction of the NS surface could be involved in the 
emission, or alternatively some fraction of the disk 
fiux could be intercepted and scattered up in a 
corona. 



16 



KS1731-260 







28/10 




29/10 




30/10 




C+L+MCD C+L+BB 


C+L+MCD C+L+BB 


C+L+MCD C+L+BB 


kTE 
r 

kTin or kTsB 
RinV COS 9 or Rbb 
EqW 


10.61°:^ 

1 Q+O-2 
-'-•O-0.2 
9 9+0.5 
^■^-0.3 


2.8tj;- 

0.9t°i 

1 0+0.1 

q + 0.5 

0.»_o.4 

136t^9 


9 7+0.1 
^•'-0.1 

12.11^7 

1 7+0.2 
4- '-0.2 
9 c+0.8 

118t30 


4-^-0.1 

q + 1.2 

d.8_Q 2 
106tp 


2 fi+'J-^ 
^•"-0.1 

1 q q+0.5 
I'-'-'-'-O.e 

o.sil^J 

1 q + 0.2 

4-3-0.1 

3.9tjj 

94li 


12.6t°J 

n 7+0.3 

^•'-0.2 

1 0+°-^ 

-■■•U-O.l 


(d.o.f) 


50.6 (51) 


48.6 (51) 


45.7 (51) 


43.9 (51) 


51.4 (51) 


51.2 (51) 


Fl-20 koV 


8.4 


8.7 


9.4 


9.6 


9.5 


9.6 


F2O-2OO keV 


0.07 


0.07 


0.08 


0.08 


0.09 


0.09 


/bb(%) 


13.5 


6.3 


13.0 


5.9 


9.8 


3.2 



Table 4: Best fit spectral results for KS1731-260 for three segments of the observation (Oct 28th, 29th, 30th). 

Nh has been set to 1.3 x 10^^ H atoms cm^^. The best fit model is the sum of a Compst model (C), a MCD 
or a BB, and a 6.4 keV gaussian line of fitted width (ape)- Rinvcos^ and Rbb are scaled at 8.8 kpc. See 
Table 1 for definition of /bb and units of Fi-20 keV and F20-200 keV- 



We have also tried the Comptt model for the 
Comptonized component. However, the tempera- 
ture for the seed photons (kT^i/) is typically 0.2- 
0.3 keV, and is therefore unreliable since the peak 
of the Wien law, at 3kTiy, is also below 2.5 keV 
(the peak is needed for kTw to be determined 
reliably). For the comptonizing cloud, kTp are 
consistent between Comptt and Compst, whereas 
T is a factor of ^ 2 lower for Comptt than for 
Compst (spherical geometry assumed). 

The broadening of the line (due to Comptoniza- 
tion?) is significant. Assuming the line is narrow 
(o-Fe=0.1 keV), we obtain a larger (87.1 for 52 
d.o.f). The hypothesis that the line is narrow is 
therefore rejected at more than 99.99% confidence 
level (j;„=40.4, P{F > F^)=5.6 x 10"^). If one 
leaves the energy of the line as a free parame- 
ter, its value tends to move towards lower energies 
(~ 6.2 keV, x^=46.1 for 50 d.o.f). However, this 
shift of the line energy is not significant (Fm=2.1, 
P{F > F™)=0.14). 

If an iron line of EqW of - 100 - 150 cV is 
present, one might expect to detect also a Comp- 
ton reflection component. First adding an edge at 
7.1 keV (cool iron) never improved the fit signif- 
icantly (significance always lower than 95%, the 
line energy and width were kept fixed). Similarly, 
we tried to fit the data with the Pexrav model, as 
a substitute to the Compst model. This never im- 
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Fig. 11. — The data and folded model (top panel), 
the residuals (center panel) and the unfolded spec- 
trum (bottom panel) of KS1731-260 for the Oc- 
tober 28th observation. The contribution from 
the Compst component is shown by a dot-dashed 
line, the blackbody is the dotted line, whereas the 
broad iron line is shown with a dashed line. 
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proved the fit either. The non detection of such a 
component might be simply related to the curvey 
shape of the spectrum (in particular the contin- 
uum is not a CPL as assumed in Pcxrav): reflec- 
tion is easier to observe when the continuum is 
PL-like (Note that if there is indeed a reflection 
component; this will affect the fitted parameters 
for the continuum and hence the strength of the 
iron line). 

The same model is also the best fit of the Oc- 
tober 29th and 30th spectra. In Table 4, we list 
the best fit parameters for the Compst+BB+Line 
and Compst+MGD +Lme models. The data and 
folded model, the residuals and unfolded spectrum 
KS1731-260 for October 28th are shown in Fig. 11. 

For KS1731-260, Li_2o kev was ~ 8.0x 10^'^ergs 

at the beginning of the observation and 
reached ~ 9.1 x 10^''ergs s~^ at the end (d=8.8 
kpc). Such X-ray luminosity is slightly larger 
than the time averaged ASM value. Throughout 
the observation, L20-200 keV contributed just at 
the level of~ 1% to the 1-200 keV luminosity of 
the source. 

4. Discussion 

We have here reported RXTE PCA+HEXTE 

timing and spectral observations of four type I 
X-ray bursters (hence NSs). In the remainder, 
first we discuss their timing properties, then their 
broad band spectral properties, and finally address 
issues concerning possible differences between BHs 
and NSs. 

4.1. Timing properties 

1E1724-3045, SLX1735-269 and GS1826-238 
were in a so-called low state (LS) whereas KS1731- 
260 was in a high state (HS). The LS sources all 
have very similar PDS, characterized by High Fre- 
quency Noise (HFN, or flat top noise. Van der 
Khs 1994) with large RMS amplitudes (30%, 2-40 
keV). 1E1724-3045 and SLX1735-269 are charac- 
terized by the presence of a ~ 1 Hz QPO, and a 
break (at i^Brcak) at 0.1-0.2 Hz. For GS1826-238, 
depending on the modeling of the PDS, fBreak ~ 
0.02 Hz or ~ 0.2 Hz, and uqpo ~ 0.2 or ~ 1.0 
Hz. KS1731-260 has a very diflFerent PDS, dis- 
playing both Very Low Frequency Noise (VLFN) 
and a QPO at ~ 7.6 Hz. These properties are all 
consistent with those of "atoll" sources at vary- 



ing luminosities; the 3 LS sources would be in the 
"island" state, whereas KS1731-260 would be in 
the "banana" state, according to the terminology 
defined by Hasinger & Van der Klis (1989). 

4.1.1. Origin of the high frequency noise? 

The high frequency noise is seen only in sources 

displaying significant hard X-ray emission (all 
sources except KS 173 1-260), and hence must be 
related to the existence of a hot scattering cloud. 
This noise in the PDS is very similar in shape and 
normalization to that seen from LS BHs, strongly 
suggesting that the same physical mechanism is 
responsible for the spectrum and variability in 
these systems. Since the spectral shape strongly 
suggests thermal Compton up scattering, then the 
variability must have something to do with vary- 
ing either the seed photons, or optical depth, or 
the dissipation in the hot region. 

The overall PDS shape can be described by 
a superposition of randomly occuring flares (or 
"shots" ), generally with some distribution of event 
durations. Such "shot noise" models can give a 
good description of the PDS (Lochner et al. 1991 
and references therein, Miyamoto et al. 1992, 
Nowak et al. 1999, Olive et al. 1998), but are 
only a phenomenological, rather than a physical, 
description of the variability. To get a physi- 
cal description requires associating the "shots" to 
changes in the comptonizing cloud. 

The simplest variability to envisage is a change 
of the soft photon input. Thermal comptonization 
involves multiple scattering of these seed photons 
on the hot electrons, so the spectrum at lower en- 
ergies responds first, followed after several scat- 
tering timescales by the higher energy spectrum. 
Such models predict that there should be time lags 
between the hard and soft energy bands, but that 
these lags (which are simply a measure of the size 
of the scattering cloud) should be constant irre- 
spective of whether the input variability was a 
short flare or a longer event. This directly conflicts 
with the observed lags in the BHCs (Miyamoto et 
al. 1988) and the NSs (Ford et al 1999), which 
clearly show longer lags for longer timescale vari- 
ability (Miyamoto et al. 1991). 

A satisfactory explanation of the long timescale 
lags is difficult. If the variability arises from vary- 
ing the seed photons then the length of the lag 
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directly implies that the region is large. This led 
Kazanas, Hua and Titarchuk (1997) to develop the 
"Extended Atmosphere Comptonization Model" 
(see also Hua, Kazanas & Cui 1999). This assumes 
a source of white noise at the center of the scatter- 
ing cloud which has a density profile n(r) oc 1/r 
out to radii of a few light seconds. Their inhomo- 
geneous density distribution appears to match the 
observed time lags, and source spectra, but the 
physical situation is very hard to envisage. The 
majority of the gravitational potential energy is 
close to the compact object, so how can this power 
a hot corona whose size is many orders of magni- 
tude larger ? Another weakness of that model, is 
that it produces different PDS at different ener- 
gies, contrary to what is observed (e.g. Nowak et 
al., 1999, Olive et al. 1998). 

An alternative model which can fit the PDS and 
time lags but with a small source size has been de- 
veloped by Poutanen & Fabian (1999). Here they 
associate the "shots" with magnetic flares above a 
cold accretion disk, and use the spectral evolution 
of the flares to produce the long lags. The flare 
begins with electrons being heated but the back- 
ground disk seed photon density is high, so the 
early time flare spectrum is soft. As the heating 
progresses, the energy dissipated in the flare dom- 
inates that of the disk, so the spectrum becomes 
hard. The drawback with such models is that they 
rely on details of the magnetic dissipation, which 
are not well known. Nonetheless, it is encouraging 
that at least under some circumstances the lags 
and PDS can be matched by the envisaged small 
source. 

4-1-2. Origin of low frequency QPOs? 

The origin of the ~ 1 Hz QPO (at vqpo) in 
the LS sources, once suspected to be a BH sig- 
nature, is also unknown, and is not accounted for 
by the above models. Chen & Taam (1994) have 
proposed that they might arise from disk lumi- 
nosity oscillations resulting from thermal viscous 
instabilities developing within the inner disk re- 
gion. Vikhlinin et al. (1994) have also suggested 
that they could result from a weak interaction be- 
tween the "shots" when the instability is triggered 
in a region of stable energy supply. To keep the 
energy released constant on large timescales, the 
appearance of a strong shot should affect the am- 
plitude and the probability of occurence of sub- 



sequent "shots". Recently, developing a model 
specific for NSs (i.e. not accounting for the sim- 
ilarities between BHs and NSs), Titarchuk & Os- 
herovich (1999) have associated the QPOs with 
radial oscillations in a boundary layer. 

It has been shown that in many sources there 
exists a strong correlation between t'Broak and 
t^QPO (Wijnands & Van der KHs 1999b). Our 
three LS sources follow this correlation (Olive & 
Barret 1999). This correlation suggests that the 
longest fluctuations (scaling as i^reak"^) and the 
QPOs are related to the same unknown physical 
mechanism or are produced in regions interacting 
with each other. Since the same correlation is ob- 
served among BHs and Z-sources, this mechanism 
cannot be related to the presence of a hard sur- 
face or even a small magnetosphere. In addition, 
in several systems (e.g. 4U1608-52, Yoshida et 
al. 1993, sec however above for SLX1735-269), a 
strict correlation between i^Break and the inferred 
mass accretion rate has been observed. For BHs, 
Esin et al. (1998) have speculated that the rapid 
variability that is associated with the hard X-ray 
emission could originate from an ADAF, with the 
variability timescales determined by the viscous 
or dynamical timescales in the ADAF. They have 
further proposed that the 1-10 Hz QPOs could re- 
sult from the interactions at the boundary (the 
so-called transition radius in ADAF terminology) 
between the ADAF and the outer cool accretion 
disk, in which case the vqpo would be some mul- 
tiples of the Keplerian frequency at the transition 
radius. Hence, a shrinking of the inner disk ra- 
dius (possibly due to an increase in the mass ac- 
cretion rate through the disk) should result in an 
increase of the QPO frequency. The correlation 
between the position of the inner disk radius in- 
ferred from the kilo- Hz QPO frequency and vqpo, 
observed for instance in 4U1 728-34 (Ford et al. 
1998), is consistent with the above picture. Fur- 
thermore, if t'Break is somehow related to the size 
of the ADAF, the observed correlation between 
'^Broak and VQpo would imply that the ADAF con- 
tracts when the accretion rate increases (possibly 
due to an increase of cooling flux from the disk). 

4-1.3. Spectral states and high frequency QPOs 

Finally, in the three sources displaying high 
frequency noise and a hard X-ray tail, no HFQ- 
POs (above 300 Hz) are detected (e.g 1E1724- 
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3045, Barrot ct al. 1999b). More generally, it 
seems that no HFQPOs are seen when i^reak is 
low (< 6 - 7 Hz; e.g. 4U1705-44, Ford et al. 
1998). The lack of HFQPOs might therefore be 
related to the presence of a hot scattering corona. 
Smearing of the HFQPO signal in such a corona 
is a possible mechanism, which is however invoked 
at higher accretion rates (i.e. for larger optical 
depth T > 5, cooler corona, Brainerd & Lamb, 
1987). This might indeed explain the disappear- 
ance of the HFQPOs after it has reached satura- 
tion (i.e. when the inner disk is at the last stable 
orbit, e.g. 4U1820-30, Bloscr ct al. 1999), or more 
generally when the source enters the upper branch 
of the "banana" state (Miller et al. 1998). The 
anticorrclation between the QPO amplitude and 
the QPO frequency observed in KS1731-260 (see 
above) is consistent with this picture (note that 
there is a weak indication that the optical depth 
derived from the spectral fitting increases when 
the source luminosity increases, see Table 4). For 
the LS sources, from the spectral analysis, we have 
derived an optical depth of a few for the Comp- 
tonizing cloud. Could that cloud, through smear- 
ing, account for the lack of HFQPOs? Using our 
upper limit, one can set a lower limit on the size of 
the scattering region. If the radiation is scattered, 
the RMS amplitude at infinity {A^c ) of a luminos- 
ity oscillation with frequency u and amplitude Aq 
at the center of a spherical region of radius Rc and 
optical depth r is A^o ~ (2^ .xe"^ + e~'^)Ao, where 
X = ^^SniyRcT/c (y is the frequency of the os- 
cillation, Kylafis & Phinney, 1989; Miller et al., 
1998). For a beaming oscillation, the attenuation 
is stronger (due to the fact that the scattering pro- 
cess tends to isotropize the photon distribution), 
and a factor 2/(1 -f r) has to be multiplied to 
the first term of the previous equation (Kylafis & 
Phinney, 1989). Let us first consider u = 300 Hz, 
from the previous expression, taking A^o = 2.5% 
(our upper limit), assuming Rq = 200 km, one 
can set an upper limit 4 and ^ 6% on for a 
luminosity and beaming oscillations, respectively. 
For a signal at 1000 Hz, these upper limits become 
12% and 20% respectively. Alternatively, if we as- 
sume = 10%, for the signal to be attenuated 
down to our upper limit, the size of the corona has 
to be larger than 100 km and 130 km for a beaming 
and luminosity 1000 Hz oscillation. Clearly, this 
indicates that with the low optical depth derived 



for the comptonizing cloud, and for any plausible 
sizes of such a cloud, the attenuation is not very 
strong, and if signal there is, it has to be intrinsi- 
cally weak. 

Another possibility which may explain the lack 
of HFQPOs in LS sources might be related to the 
position of the inner disk radius. Obviously, if 
the HFQPOs are generated in the disk, and if the 
inner edge of the disk lies at large radius (see be- 
low), no signals at such high frequencies should 
be seen. In that respect, it is worth pointing 
out that the bump (around 10 Hz) in the PDS 
of 1E1724 3045 (Fig.6) was interpreted by Psaltis 
et al. (1999) as a low frequency kilo-Hz QPO (the 
same bump is seen in GS1826-238 and SLX1735- 
269). Their conclusion was derived from the fact 
that the 0.8 Hz QPO and the bump at 10 Hz fitted 
in the global correlation observed over several fre- 
quency decades between i^qpo and the frequency 
of the lower kilo-Hz QPOs (Psaltis et al. 1999). 

4.2. Spectral properties 

To first order, the broad band LS spectra of 
NSs can be approximated with a PL followed by 
an exponential cutoff at 50-80 keV. There is also 
a weak soft component in 2 of the 3 LS sources, 
with temperature of ~ 0.5 — 1 keV contributing 
less than ~ 30% of the total luminosity. The lat- 
ter quantities, when derived from the fit of PCA 
data are subject to uncertainties because the fit 
starts at 2.5 keV, because they depend on the 
Nfj values assumed, and because there are some 
calibrations uncertainties of the PCA at low en- 
ergies (E < 5 keV). However, as for example in 
the case of 1E1724--3045, the values derived from 
our spectral fitting are consistent with the ones 
obtained from observations that are more sensi- 
tive to such a soft component (e.g. by Beppo- 
SAX or ASCA, Guainazzi et al. 1998, Barret et 
al., 1999a). Thermal Comptonization of these soft 
photons by hot electrons then appears to be the 
most plausible emission mechanism for these sys- 
tems, and for a spherical scattering cloud, the de- 
rived optical depths and electron temperatures are 
in the ranges of 2 < t < 4 and 15 < kTg < 30 keV. 
Thermal Comptonization also dominates the spec- 
tral formation in KS1731-260, but with a signifi- 
cantly lower kTE (~ 3 keV) and larger r (~ 10). 
The main difference between a source in the HS 
and a source in the LS is best illustrated in Fig. 
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Fig. 12.— A uFu plot of the spectra of KS1731- 
260 (high state) and GS1826-238 (low state). For 
KS1731-260, the bulk of the energy is radiated be- 
low ~ 7 keV. For GS1826-238, on the other hand, 
its spectrum is relatively flat, peaks around 40 keV 
and drops rapidly above 100 keV. No scaling for 
the different source distances has been performed. 

12, where a uFu plot of the PCA/HEXTE spec- 
tra of KS1731-260 and GS1826-238 are shown. 
Clearly for KS1731-260, the bulk of the energy 
is radiated below 10 keV, whereas for GS1826- 
238 the energy spectrum is flatter and a large frac- 
tion (~ 50%) of the energy is radiated in the hard 
X-ray band. 

4-2.1. Iron line emission and reflection in NS 
LMXBs 

There is strong evidence for iron line emission 
in 2 out of the 3 LS sources examined here (the 
evidence is weaker for 1E1724 3045). The derived 
line energy is consistent with 6.4-6.5 keV, i.e. flu- 
orescence from neutral/moderately ionised iron. 
There are three possible sites for line emission in 
LMXBs, the NS surface itself, the accretion disk, 
and accretion disk coronal wind. Any line from the 
NS would be strongly redshifted, down to ^ 5 kcV, 
while the accretion disk corona should be strongly 
ionised, producing iron lines at 6.7 and 7.0 keV. 
Thus, the most likely origin for the observed 6.4 
keV line is irradiation of the accretion disk by the 
X-ray source. If so, then it should be accompa- 
nied by a Compton reflected spectrum. We sig- 



nificantly detect a neutral reflected continuum in 
GS1826-238, and a neutral/moderately ionised re- 
flected continuum in SLX1735-269. 

The reflection albedo, the ratio between the re- 
flected and incident luminosities, is then an ob- 
servationally determined quantity from our data. 
The reflection probability is determined by the rel- 
ative importance of electron scattering and photo- 
electric absorption. At low energies the X-ray 
photons tend to be photo electrically absorbed 
rather than scattered, so the reflection probability 
is low, while at higher energies the opposite is true. 
The photo-electric absorption probabilities are de- 
termined by both the element abundances (more 
heavy elements mean more hard X-ray opacity 
and so less reflection) and ionisation state (high 
ionisation states mean less bound electrons, so less 
opacity and more reflection). Beyond 10-20 keV, 
the photo electric opacities become negligible in 
comparison to the electron scattering probability, 
but at these energies the electron scattering is not 
elastic. The incoming photon can lose much of its 
energy to the electron as it scatters, leading to a 
marked decrease in the reflected probability above 
50-100 keV. Thus, even high energy photons can 
dump most of their energy in the disk. For hard 
spectra, such as those seen in GS1826-238, the 
maximum albedo, even with complete ionisation of 
the disk is ~ 0.75, while the observed albedo is 
much less, at ~ 0.25, because the ionisation state 
of the reflecting material is low/moderate, so most 
of the low energy incident flux is photo-electrically 
absorbed rather than reflected. For the HS spec- 
tra, such as KS1731-260, the reflection albedo can 
be much higher if the disk is completely ionised, 
since the spectrum is much softer. This is im- 
portant because the disk heating and hence opti- 
cal/UV reprocessed X-ray flux is determined by 
the non reflected flux, {1 — ax). Various attempts 
to derive the albedo from the optical reprocessed 
flux (e.g. de De Jong, van Paradijs & Augusteijn 
1996) give ~ 0.9. This is not in conflict with 
our values since most of the LMXBs considered 
were in the high state, where much more of the 
spectrum can be reflected. However, we caution 
that the optical determinations assume that the 
disk shape is that given by Vrtilek et al. (1990), 
with height oc r^^"^. This was derived assuming 
that the disk is isothermal with height (central 
temperature equal to surface temperature) , which 
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is incorrect (e.g. Dubus et al. 1999). 

While the fraction of incident to reflected flux 
depends on the X-ray albedo (which in turn is 
given by the elemental abundances and ionisation 
state of the reflector and the spectral shape of the 
incident spectrum), the total amount of reflection 
seen depends on the; solid angle subtended by the 
material, 0/27r =frefl, and its inclination to the 
line of sight. For an assumed inclination angle of 
60°, we derived frcfl= 0.15 and 0.28 in GS1826- 
238 and SLX1735-269, respectively, with an asso- 
ciated iron line EqW of 50 and 39 eV. The line 
EqW and especially the reflected fraction arc sig- 
nificantly lower than the ~ 130 eV and frefl= 1 
for an isotropic X-ray source above a flat, inflnite 
disk (George & Fabian 1991; Matt, Perola & Piro 
1991). 

4-2.2. Accretion Geometry 

Both the broad band continuum spectrum, and 
its variability (e.g. van der Klis 1994; Ford et al. 
1999) and the properties of the reflected spectrum 
(Zycki, Done & Smith 1998, 1999) arc very similar 
to that seen in the LS BH systems, suggesting that 
the same physical mechanisms operate in a simi- 
lar accretion gcomc^try. If so. then the hard emis- 
sion cannot have anything to do with the magne- 
tosphere or the NS surface, and must instead be 
connected with the accretion flow. Mechanisms 
suggested for the hard power law in BH systems 
which can also work in NSs are either an accre- 
tion disk corona (magnetic reconnection in flares 
above a disk), or a hot accretion flow. The lack of a 
strong reflection signature and weak soft emission 
is at first sight incompatible with the magnetic 
flaring model. A geometry where the hard X-ray 
source is above an optically thick disk should lead 
to ffofi ^ 1 rather than the frofi^ 0.3 observed here. 
In addition, of the hard X-ray fiux, half will go 
up and escape while half will go down and illumi- 
nate the cool material. Of these, 10% can be re- 
flected as a hard X-ray component (if the material 
is mostly neutral, and has solar abundances), but 
the remaining 90% arc thermalised, and shoifld 
emerge as a soft component. The soft photon lu- 
minosity would then be roughly half of that of the 
hard component. However, it is only emitted into 
27r solid angle because it is optically thick, and 
since the hard component is emitted into 47r; the 
inferred luminosities of the soft and hard compo- 



nents should be about equal (see Gierlinski et al. 
1997). 

The properties of the reflected component and 
line emission together with the observation of time 
smearing or delays in the optical counterparts of 
X-ray bursts testify of the presence of an outer 
disk in those systems. It would seem however that 
a significant fraction of the disk cannot be seen, 
either in reflection, or by reprocessing the hard 
X-ray radiation, or by emission. So the question 
is: is there then an inner disk at all, say within 
~ 50 — 100 Rs ? If not, then the disk corona 
model can be ruled out. However, the disk reflec- 
tion and reprocessing signatures could perhaps be 
masked by a very highly ionised inner disk (Ross 
et al. 1998, but see Done & Zycki 1998), while 
the intrinsic disk emission is then cooler and less 
intense if the power is mostly dissipated in the 
corona (Svcnsson & Zdziarski 1994). Another pos- 
sibility is that the hot plasma has a bulk veloc- 
ity f3c directed away from the disk (Beloborodov 
1999). The bulk motion might be due to the pres- 
sure of the reflected radiation and/or plasma ejec- 
tion from magnetic flares. A mildly relativistic es- 
caping flow (/? ~ 0.3) causes aberration which re- 
duces the downward irradiating flux which in turn 
reduces the feedback of reflection and reprocess- 
ing. 

Perhaps the most telling discrepancy in all 
these models is that there should also be strong 
emission from a boundary layer where the disk 
and NS interact. The boundary layer emission can 
be twice as large as that from the accretion disk 
(whether standard disk or disk powered corona) 
once relativistic corrections arc included (Shakura 
& Sunyaev 1986). There is no strong soft emis- 
sion component in the spectrum, ruling out an 
optically thick boundary layer (unless it also dissi- 
pates most of the accretion energy in an optically 
thin region perhaps powered by magnetic recon- 
nection). King & Lasota (1987) calculated that 
the boundary layer between a standard disk and 
the NS can be optically thin only for mass accre- 
tion rates much lower than those considered here 
(but see recent calculations by Inogamov & Sun- 
yaev, 1999 and Popham & Sunyaev 1999). Alter- 
natively, if the disk truncates at the last stable or- 
bit (allowed only for certain NS equations of state) 
then the boundary layer is between material freely 
spiraling in rather than a disk. This can give an 
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optically thin boundary layer up to luminosities 
of a few percent of Eddington, as required here 
(Kluzniak & Wilson 1991). However, it seems un- 
likely that an optically thin, hot boundary layer 
could be similar in temperature and variability to 
disk emission produced by a completely different 
process of magnetic flare. Removing the boundary 
layer by ejecting the material just before it reaches 
the NS surface (e.g. by a centrifugal magnetic field 
barrier or winds) is not a viable way around these 
difficulties. The observation of X-ray bursts from 
these systems proves that the accreting material 
does accumulate on the NS surface (e.g. GS1826- 
238). 

Such a series of requirements make the disk 
corona models seem rather contrived, so we in- 
vestigate the alternative model where the inner 
disk is replaced by an X-ray emitting hot flow. 
This would naturally account for the weakness of 
reflection and reprocessing, and the relative weak- 
ness of the intrinsic soft emission. It also provides 
an explanation for the strong correlation observed 
between the relative strength of the reflection com- 
ponent (frofl), and the photon index of the intrin- 
sic PL (r) in Seyfert galaxies and galactic BHs 
(Zdziarski et al. 1999). The same correlation also 
applies to NSs, including those considered here; 
for GS1826-238, frcfl=0. 151^4 and r=1.72t°;^^, 
whereas for SLX1735-269, the corresponding val- 
ues are frcfl=0.28ti^:'i'^ and r=2.09t2 Jj^). This cor- 
relation can be understood by considering a hot X- 
ray plasma that fills a variable-radius hole in the 
inner disk; as the disk radius decreases, its solid 
angle increases, thereby increasing both the reflec- 
tion and the flux of cool cooling photons. The 
stronger cooling of the hot plasma then leads to a 
steepening of the PL slope (Zycki, Done & Smith 
1998, Zdziarski ct al. 1999)^. 

The above scenario, where the innermost radius 

'Alternatively, this correlation can also be explained by the 
outflow model of Beloborodov (1999). Increasing the semi- 
relativistic outflow velocity, /9, leads to a decrease of both 
r and frefl, and yields a correlation that seems to fit better 
the observations than the two phase disk model (Zdziarski 
et al. 1999). Another interesting feature of the ejection 
model, is that it can explain i^^d > 1 (as observed in some 
Seyferts) because for negative values of /?, the coupling be- 
tween the ejected plasma and the disk can be very strong. 
As said above, however, the weak point of such a model is 
that it relies on the poorly known physics of the magnetic 
dissipation within the disk. 



of the accretion disk moves as a function of spec- 
tral state, can also naturally explain the strong 
correlation between spectral properties and the 
frequency of the HFQPOs (e.g. Kaaret ct al. 1998; 
Mendez et al. 1999, Bloser et al. 1999). While 
there are several current models proposed for these 
high frequency features in the power spectrum, 
they all require that the NS spin beats with an in- 
ner disk frequency (see e.g. the review by Van der 
Klis 1998). The easiest way to change the inner 
disk frequency is to change the inner disk radius. 
This then naturally explains the lack of HFQPOs 
in the LS NSs as being systems where the inner 
disk is truncated at large radii, where the inter- 
action with the NS surface is small. As the disk 
moves inwards, the spectrum softens (F increases; 
e.g. in 4U0614+09 and 4U1608-52, Kaaret et 
al. 1998, and the disk contribution increases; e.g. 
4U0614+09, Ford ct al., 1997), the HFQPOs ap- 
pear, and their frequencies increase (Piraino et al., 
1999). 

How can we reconcile our spectral and timing 
observations in that picture? GS1826-238 would 
have an accretion disk with a large inner radius. 
Its inner disk temperature is then too low for 
its contribution to be detected with the PGA. 
The cooling is low, its spectrum is hard and the 
amount of reflection is rather small. Conversely, 
for SLX1735-269, the disk extends closer to the 
irradiating source, the solid angle for reflection 
is larger and its spectrum is softer. Finally for 
KS1731-260, the disk gets even closer; kilo-Hz 
QPOs are produced and are seen before disappear- 
ing as a response to a further increase of the ac- 
cretion rate. The cooling is strong and leads to a 
quenching of the hard X-ray emission; at the same 
time the iron line gets stronger (so probably, does 
the reflection component). The above picture, al- 
though attractive, does not explain why 1E1724- 
3045, which has spectral and timing parameters 
more or less similar to SLXl 735-269 does not dis- 
play any reflection component (yet it shows a weak 
iron line) . One possible explanation could be that 
the inclination is larger for 1E1724-3045 than for 
SLX1735-269. 

4-. 2. 3. Emission Mechanisms 

But what can replace the inner accretion disk ? 
The recent rediscovery of a stable. X-ray hot so- 
lution of the accretion flow equations has caused 
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much excitement. The main assumption of these 
advective solutions is that the gravitational en- 
ergy released by viscosity is gained mainly by 
the protons, and that these heat the electrons 
only by Coulomb collisions. For low mass accre- 
tion rates the flow is optically thin so that the 
Coulomb collision rate is very low. The proton 
temperature is high, so the flow has a large scale 
height (quasi-spherical). The small amount of en- 
ergy that is gained by the electrons is radiated as 
cyclo/synchrotron emission, bremsstrahlung and 
Comptonization of these seed photon distributions 
(Narayan & Yi 1995). Pure advective models have 
a disk existing only at very large radii, where its 
emission is negligible. The resulting spectrum is 
then rather hard, since the only seed photons are 
the self-produced ones in the optically thin flow, 
and the electron temperature is of order 100 
keV (Narayan & Yi 1995). As the mass accre- 
tion rate increases, the flow becomes denser, so 
the Coulomb collisions are more effective at trans- 
fering energy from the protons to the electrons, so 
the radiative efficiency increases. But this process 
cannot continue indefinitely: as the flow becomes 
optically thick, the Coulomb collisions drain all 
the energy from the protons, and the flow collapses 
into a standard, optically thick, geometrically thin 
accretion disk (Esin et al. 1998). 

However, there is a clear difference between ad- 
vective models for BH and NS systems. For BHs 
the energy advected with the protons in the flow 
can be swept invisibly down into the black hole. 
Such a flow will collapse at i 0.4a^ LEdd where 
OL ^ 0.2 (Yi et al. 1996) is the a disk viscos- 
ity. For NSs the advected energy is released in 
a boundary layer as the flow hits the NS surface. 
If this boundary layer is optically thick, then the 
increase in seed photons for the Compton cooling 
will cause the advective flow to collapse at Edding- 
ton scaled mass accretion rates lower (by at least a 
factor 3) than for BHs (Narayan & Yi 1995; Yi et 
al. 1996). Yet the spectral state transition takes 
place at roughly the same Eddington fraction for 
NSs and BHs, namely at ~ 0.05 — 0.1 Lsdd (see 
e.g. Mitsuda et al. 1989 for the NS 4U1608-52, 
and for the BH Cyg X-1, Esin et al. 1998, see also 
Esin et al. 1997), and no strong soft emission from 
a boundary layer is seen. 

One way out of this impasse is if the bound- 
ary layer is optically thin. Explicit calculations 



for the boundary layer between an ADAF and NS 
surface have not yet been done, though rather dif- 
ferent spherical flows have been shown to give op- 
tically thin boundary layer emission up to fairly 
high luminosities (e.g. Zane et al. 1998). Thus 
an optically thin boundary layer between an ad- 
vective flow and NS seems plausible. In this case, 
the boundary layer is an additional heating source 
for the hot plasma in the advective flow, so it 
might be expected to merge rather smoothly in 
both spectral and variability properties with the 
rest of the hot, optically thin accretion flow. The 
extra photons from the boundary layer are then 
hard. This lowers the temperature of the advec- 
tive flow slightly, but the critical mass accretion 
rate which can be sustained in the flow is then 
very similar to that calculated for the BH case 
(Narayan & Yi 1995). The one caveat to this 
is that even if the boundary layer were optically 
thin then the NS surface would intercept and ther- 
malize some of the boundary layer emission, again 
leading to an additional source of soft seed pho- 
tons. If the NS surface is ionised and/or mainly 
made up of hydrogen due to settling of heavy ele- 
ments, then the reflection albedo can be as high as 
0.6. Thus at least 25% of the boundary layer lumi- 
nosity should emerge as a soft component. We are 
currently calculating the effect of such an optically 
thin boundary layer plus its thermalized emission 
on the properties of the advective accretion flow, 
to see whether it can allow the observed similarity 
in state transition accretion rate between BHs and 
NSs (Done & Barret 1999). 

It has been recently suggested that mass loss 
via winds was a natural consequence of ADAFs, 
and was such, that only a tiny fraction (^ 1) of 
the gas siipplied through the ADAFs was ac;tually 
accreted onto the central object (Blandford and 
Begelman 1999, Quataert & Narayan 1999). Such 
winds would provide an alternative explanation 
for the dimness of quiescent BHs, attributed in 
the ADAF models, to low radiative efficiency and 
the existence of event horizons in those systems 
(Narayan, Garcia & McClintock 1997; Menou et 
al. 1999). However, in our picture, where the NS is 
surrounded by an ADAF, the observation of X-ray 
bursts tells us that all the matter flowing through 
the ADAF can be accreted onto the NS. So, un- 
less the mass transfer rate is very much larger 
than currently thought in LMXBs (see Menou et 
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al. 1999 for a recent discussion), the presence of 
ADAFs around bursting NSs would strongly ar- 
gue against the existence of powerful winds from 
such accretion flows (i.e. ADIOS in Blandford & 
Begelman 1999). 

To summarize, the similarity of the BH and NS 
spcictra and variability in their LS strongly suggest 
that the same mechanisms are operating i.e. that 
the hard X-ray emission is connected to the accre- 
tion flow witli the NS surface having little impact 
on it. Although the latter remains an unescapable 
part of the system, there is nothing which can 
be firmly associated witli the strong emission ex- 
pected from the surface or the boundary layer even 
though it is known that the accretion material has 
to accumulate onto the NS surface because of X- 
ray bursts. This problem applies to both the disk 
corona and advective flow models. The bound- 
ary layer emission then must either be optically 
thin, or emitted at too low a temperature to be 
observed. The latter seems unlikely given pre- 
vious ASCA and Beppo-SAX observations, hence 
we conclude that the boundary layer is optically 
thin, and hot (see recent calculations by Popham 
& Sunyaev, 1999, Inogamov & Sunyaev, 1999). 
However, it seems difficult to imagine a situation 
where the boundary layer emission could have sim- 
ilar spectral and variability properties to magnetic 
flares above a disk. Hence we favor models where 
the inner disk is replaced by an X-ray hot, opti- 
cally thin flow. The only known stable hot solu- 
tions to the accretion equations are the advective 
flows. Unfortunately ADAF solutions have not 
been developed with self-consistency for the NS 
case, and it is not yet known whether these can 
have an optically thin boundary layer. If so then 
this emission could well merge smoothly with the 
hot optically thin advective flow. We stress that 
the maximum mass accretion rate at which an ad- 
vective flow can be sustained holds out the pos- 
sibility of showing whether or not advective flows 
can really be present. Observations show that the 
hard/soft spectral transition associated in these 
models with the collapse of the advective flow oc- 
curs at roughly the same mass accretion rate in 
both NS and BH systems. While an optically 
thick boundary layer produces a large difference 
in this critical mass accretion rate, an optically 
thin boundary layer does not. However, even an 
optically thin boundary layer produces some soft 



photon flux from reprocessing of the flux illumi- 
nating the NS surface, which will lower this crit- 
ical mass accretion rate. Further calculations are 
needed to see whether this effect is small enough 
to match the observed spectral transitions in NSs. 

4.3. Comparison between BHs and NSs 

Thanks to RXTE and Beppo-SAX, the number 
of LS NSs observed simultaneously in X-rays and 
hard X-rays is growing rapidly, so that reliable 
comparisons between BHs and NSs can now be 
carried out. Although, as illustrated in this paper, 
the most recent data indicate that BHs and NSs 
are very similar in many respects, especially in 
their low states, it remains critical to search for 
observational criteria that could distinguish these 
two types of accreting systems. 

4-3.1. Spectral shape differences? 

Heindl & Smith (1998) have pointed out that 

the index of the power law part of the spectrum is 
signiflcantly larger for NSs than for BHCs (r> 1.8 
for NSs versus 1.4 <r< 1.6 for BHs, see however 
Churazov et al. 1997 for an opposite conclusion). 
The data presented here are in general terms con- 
sistent with this claim, but indicates that the sepa- 
ration is by no means large (the index for GS1826- 
238 is 1.7). On the other hand, Heindl & Smith's 
claim, is inconsistent with several recent observa- 
tions of NSs, as for example those of S AX J 1748. 9- 
2021 in NGC 6440 (In't Zand et al. 1999) or those 
of the two dippers 4U1915-05 (Church et al. 1998) 
and XB1323-619 (Balucinska-Church et al., 1999). 
In the flrst case, a fit by a broken PL yields a pho- 
ton index of 1.54 ±0.03 and 2.13 ±0.04 below and 
above the break at 18.1 ±1.2 keV (In't Zand et 
al. 1999). For 4U1915-05, the broad band 0.2-200 
keV non-dip Beppo-SAX spectrum can be fitted by 
a blackbody plus a cutoff PL of index 1.6 ± 0.01 
(Ecutoff= 80.4 ± 10 keV). Similarly, for XB1323- 
619; Balucinska-Church et al. (1999) have found 
that the non-dip spectrum is a cutoff power law 
with r= 1.48 ± 0.01 and Ec„toii =44.1 keV. Based 
on these observations, we therefore conclude that 
the above criterion is not valid. 

It has also been proposed that, in the frame- 
work of thermal Comptonization models, the elec- 
tron temperature of the scattering cloud (kTs) ap- 
peared to be systematically lower for NSs than for 
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BHs; kTE< 30 keV versus kTE> 50 kcV (Tavani 
&: Barret 1997, Zdziarski et al., 1998, Churazov et 
al. 1997). This naturally reflects the standard pic- 
ture that on average BH spectra arc harder than 
NS spectra, and has been tentatively explained by 
the additional cooling provided by the NS surface, 
which may act as a thermostat capable of limit- 
ing the maximum kTs achievable in these systems 
(Kluzniak 1993, Sunyaev and Titarchuk 1989). 
The data presented in this paper arc certainly con- 
sistent with this criterion. However, it is worth 
noting that there are some speculative BHCs (e.g. 
GRS1758-258) for which kTg derived from the fit- 
ting of their hard X-ray spectra with the Sun- 
yaev & Titarchuk (1980) Comptonization model 
{Compst in XSPEC) is below 50 keV (kTE-33 keV 
in Mandrou et al. 1994 from a fit of SIGMA data 
alone). However, as pointed out by Zdziarski et al. 
(1998), the latter temperature is probably a gross 
underestimate, and should be reevaluated using 
broad band spectra and more appropriate rela- 
tivistic Comptonization models including Comp- 
ton reflection and relativistic effects (e.g. Pouta- 
nen & Svensson 1996). If such an underestimate 
has indeed been observed in Cyg X-1 (kTE in- 
creased from 27 keV with Compst up to ~ 100 
keV with CompPS; Gierlinski et al. 1997), we note 
that in the case of 1E1724 3045, kTE derived with 
Comptt is 30 keV, whereas it is only slightly larger 
(^ 35 keV) with the relativistic CompPS model 
(the same is true for GS 1826 238). Bearing this 
in mind, it remains that this criterion should be 
considered seriously for sources showing evidence 
for thermal cutoffs in their hard spectra. More 
data should tell us soon whether all BHCs fltted 
with relativistic Comptonization models will in- 
deed have kTE ^ 50 keV, whereas all NSs fitted 
with the same models will have kTs below the 
above value. 

In any case, the above criterion applies when 
an energy cutoff is observed. However such cutoffs 
are not always present in NS hard X-ray spectra. 
The first SIGMA observation of 1E1724 3045 re- 
vealed a non attenuated hard power law (r=1.8) 
extending up to 200 keV (Barret et al. 1991). Aql 
X-1 was also observed by BATSE from 20 kcV 
up to ~ 100 keV with F in the range 2.1-2.6 and 
no evidence for a high energy cutoff (Harmon et 
al. 1996). Finally, a recent Beppo-SAX observa- 
tion of 4U0614-I-09 has placed a lower bound on 



any exponential cutoff of the power law (F=2.3) 
of 200 keV^ (Piraino et al. 1999). This means 
that just as there are two classes of BHs based 
on their hard X-ray spectra (Grove et al. 1998), 
there might also be two classes of NSs. Members 
of the first class would display hard X-ray spec- 
tra with energy cutoffs, which would result from 
thermal Comptonization. Members of the second 
class would have non-attenuated power laws (up 
to an energy which remains to be accurately de- 
termined), similar to the power laws observed in 
the soft state of BHs. Such power laws could 
be produced by non-thermal Comptonization, i.e. 
Comptonization on non-thermal particles (for re- 
cent reviews see Poutanen 1999, and Coppi 1999). 
Alternatively, Titarchuk et al. (1997) have pro- 
posed that these power laws could be produced 
by bulk-comptonization in a convergent accretion 
flow. The same mechanism seems to be railed out 
for NSs because it produces spectra much harder 
than that observed (r<l, Titarchuk et al. 1996). 
BHs with power laws have typical F in the range 
2.5-3.0, or similar to NSs in the hard X-ray band 
(e.g. 1E1724-3045). Therefore, the indication of 
a steep hard X-ray (E>30 keV) spectrum with 
F> 2.5 cannot alone be used to claim that an un- 
known system contains a NS. Similarly, a hard X- 
ray power law spectrum with r< 2.5 is not unique 
to BHs (e.g. 4U0614-F091). 

4-3.2. Luminosity differences? 

The idea that BHs and NSs are actually hardly 
distinguishable by their broad band spectral shape 
has led Barret et al. (1996) to propose a luminos- 
ity criterion (see also Barret & Vedrenne 1994, Van 
der Klis & Van Paradijs 1994). They compared 
the 1-20 keV luminosity (Li_2o kev) to the 20-200 
keV luminosity (L20-200 kov) for all secure BHBs 
(by secure we mean BHs with mass function esti- 
mates indicating a compact object of mass larger 

^Alternatively, the Beppo-SAX spectrum could be modeled 
as thermal Comptonization {CompPS model in XSPEC, 
Poutanen & Svensson, 1996) for which a very high kTg was 
derived (246^3° kcV, Piraino ct al. 1999). Such a value, 
which puts the cutoff energy above 700 kcV, is itself well 
above the high energy threshold of the Beppo-SAX/Y'DS 
(the last significant data point of the spectrum is at 180 
kcV), and therefore must be considered with cautions. 
However, if this result is confirmed by observations per- 
formed at even higher energies, it would maJje the kTg cri- 
terion discussed above definitely unvalid. 
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than 3Mq, hereafter BHBs for BH Binaries), and 
all NSs detected up to at least 100 keV. Figure 
13 is an updated version of Fig. 1 in Barret et 
al. (1996). There arc now 16 NSs detected at 
100 keV with (quasi-)simultaneous coverage in X- 
rays; all but Cen X-4 detections came within the 
last 8 years. It includes 6 more NSs recently de- 
tected by either RXTE (PCA+HEXTE) or Beppo- 
SAX (MECS+LECS+HPGSPC+PDS); namely 
the milhsecond pulsar SAXJ18G8.4-3658 (Gilfanov 
et al. 1998, Heindl & Smith 1998), the dippers 
4U1915-05 (Church et al. 1998) and XB1323-619 
(Balucinska-Church et al., 1999), SAXJ1748.8- 
2021 (In't Zand et al. 1999), SLX1735-269 and 
GS1826-238 (this work). For the two dippers, 
Li-20 keV and L20-200 kcV have been computed 
from their non-dip spectra. In addition, the 
luminosities for 1E1724-3045 (this work), for 
4U0614+G91 [Beppo-SAX observation, Piraino et 
al. 1999), and Aql X-1 (joint nearly simultane- 
ous BATSE & ASCA observations, Rubin et al. 
1999) have been updated. The figure also includes 
GRS1009-45 (Nova Velorum 1993), recently shown 
to be another secure BHB (Filippenko et al. 1999). 

For both SAXJ1808.4-3658 and 4U1915-05, we 
have assumed the distance inferred from X-ray 
burst studies: 4 and 9.3 kpc respectively (In't 
Zand et al. 1998, Yoshida 1992), and for XB1323- 
619 a distance of 10 kpc (Balucinska-Church et al., 
1999). SAXJ1748.8-2021 is located in the globu- 
lar cluster NGC6440 whose distance is estimated 
to be 8.5 ± 0.5 kpc (Ortolani et al. 1994). For 
GRS1009-45 we have estimated the distance ac- 
cording to the method described in Barret et al. 
(1996). We determine the radius of the secondary 
assuming that it fills its Roche lobe (the radius is 
then given by the orbital period Porb=6.86 hours, 
and its mass, which we take to be 0.5 M©). For 
a K7 secondary, from the absolute visual flux, we 
determine the absolute visual magnitude Af„ ^ 
8.2, according to the values tabulated by Popper 
(1980). Finally, using the apparent dereddened 
magnitude corrc;c;te<:l for interstellar reddening and 
for the contribution (fdisk) by the accretion disk 
to the continuum flirx [Vquies = 21.4 - 21.9, E(B- 
V)=G.2, fdisk=6G%, Shahbaz & Kuulkcrs 1998), 
we obtain a distance in the range 5.0 to 6.5 kpc. 
Shahbaz and Kuulkers (1998) using an empirical 
linear relationship between the orbital period, the 
optical outburst amplitude magnitude derived a 
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Fig. 13. The hard X-ray (20-200 keV) versus 
the X-ray (1-20 keV) luminosities of Black Hole 
Binaries (filled symbols) and Neutron Star Bina- 
ries (open symbols). Only BHs with mass func- 
tions indicating a mass for the compact object in 
excess of ^Mq are considered. For NSs, the lumi- 
nosities are computed from observations in which 
the source was detected up to at least 100 kcV. For 
KS1731-26G, the luminosities drawn correspond to 
the SIGMA detection (Barret et al. 1992). This 
is an updated version of Fig. 1 in Barret et al. 
(1996) from which most luminosities can be re- 
trieved. The so-called X-ray burster box is plotted 
as a dot-dashed line. Its boundaries are simply 
defined as to contain all the NSs. 

distance of 4 kpc. We adopt in Fig. 13 a distance 
of 5 kpc. The X-ray and hard X-ray fluxes for 
GRS 1009-45 have been taken from Barret et al. 
(1996) from observations reported by Sunyaev et 
al. (1994). 

Looking at Fig. 13, one can see that the dis- 
tinction in luminosity clearly holds: all NSs lie in 
the so-called X-ray burster box, whereas all BHs 
are found outside. This figure shows that compar- 
isons of L20-200 kcV and Li_20 koV can be used to 
distinguish between BHs and NSs. First, although 
the presence of a hard tail is not a unique signature 
of BH accretion, the high luminosity in the hard 
tail appears to be such a signature: only BHBs are 
observed with L2o_2oo kev^ Lcrit, where we define 
Lcrit ^ 1.5 X IG^''' ergs s~-^. Second, the hard X- 
ray luminosity of NSs can approach L^rit , but only 
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when Li_2o keV^ Lcrit- On the other hand BHBs 
can have both Li_2o kev and L20-200 kev above 
Lcrit (although L20-200 keV may then be substan- 
tiaUy lower than Li_20 kov)- Finally, the most re- 
cent broad band measurements of LS NSs as those 
presented in this paper all indicate that these sys- 
tems have a ratio Li_20 kov/L20-200 koV iu a nar- 
row range of ~ 1.0 — 2.0. This yields the corre- 
lation between Li_2o keV and L20-200 keV seen in 
Fig. 13. This suggests that below a luminosity 
~ 2x Lcrit, the spectral shape of a given NS does 
not depend much on its luminosity. However, ob- 
servation of softer spectra (almost blackbody-like, 
Asai et al. 1998) from quiescent NS transients im- 
ply that this correlation breaks down somewhere 
between ^ lO^** and 10'^^ ergs s"^. Another inter- 
esting feature of Fig. 13 is the fact that the BHBs 
(those lacking an ultra soft component) lie roughly 
on an extrapolation of the sequence formed by the 
NSs points. This again points to the underlying 
similarity of their emission sources and basic ac- 
cretion geometries. 

5. Summary 

We have reported RXTE timing and spectral 
observations of four type I X-ray bursters. The re- 
sults presented here nicely illustrate that beside its 
unprecedented timing capabilities, RXTE has also 
sufficient sensitivity for detailed spectral studies. 
Our spectral fitting has revealed that the broad 
band spectra of low state NSs are often more com- 
plicated than previously thought. In particular, 
we have detected a Compton reflected component 
plus an iron Ka line in two systems. In both cases, 
the data indicate that the most likely site for the 
reflector is a cool accretion disk, which is trun- 
cated somewhere to offer a relatively small solid 
angle to the irradiating source, which could be a 
quasi-spherical hot advection dominated inner ac- 
cretion flow. Similar behavior is observed among 
BHs leading to very similar conclusions. In addi- 
tion, within the limitations of the PCA calibration 
uncertainties at low energies, we have detected a 
soft component, which is equally well described 
by a multi-color disk blackbody or a blackbody. 
In our picture, this soft component most probably 
originates from the tnmcatcd accretion disk. 

At higher energies, the HEXTE sensitivity is 
sufficient to accurately locate the energy cutoff in 



the low state hard X-ray spectra of these systems. 
When combined with PCA data, we have shown 
that the spectral steepness of the hard tails ob- 
served by SIGMA/BATSE was not intrinsic, but 
was rather due to the presence of a high energy 
cutoff (~ 50 — 80 keV). The cutoff power law spec- 
tral shape observed strongly suggests that thermal 
Comptonization is the dominant emission mecha- 
nism in these systems. The Comptonization pro- 
cess would take place in an optically thin hot 
boundary layer merged with the central advective 
corona, located between the truncated accretion 
disk and the NS surface. 

Of the most recent criteria that have been pro- 
posed to distinguish BHs from NSs (out of quies- 
cence), two appear consistent with all the avail- 
able data. The first one states that for NSs dis- 
playing thermal Comptonization, they seem to be 
unable to achieve electron temperature as high as 
BHs, and hence Comptonization with kTE > 50 
keV would be a BH signature. This criterion is 
weakened however by the fact that there are some 
NSs that do no display high energy cutoffs in their 
hard X-ray spectra. The second one is a luminos- 
ity based argument and claims that only BHs are 
capable of emitting bright hard X-ray tails with 
luminosity larger than ~ 1.5 x 10^'' ergs s~^. 

The RXTE archive contain a wealth of observa- 
tions of NSs and BHs in different luminosity states, 
and represent a very uniform data set, that should 
be used to test the picture proposed in this paper. 
It is also clear that the same data can be used 
to challenge the observational criteria discussed 
above to distinguish BHs from NSs. 
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